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1 

l-^*?©fcfe©|| 1 ©x-<i' hJUT^^-if^iDm-r'S.;^ 

7 hS!|l©fci*©»2©;^-^i' hJU-r-5f^iR*-rs;^ 
-eti-en©SfR3n/c-9-:/TH'K*jW 

mam i ©^H«:*f L/-cffs®ai;?3nfc^rai(cffis© 
mm. 2 ] mmi>:^T v zfifi^ m(Dv^7 u -f 

i. ffit©1>-:/TU-ffiia{C43W*:*-7-fey h4S«66 

m 1 K:iet£©X'>ii' hJl'T^" ^^^(Ofc^ox^m. 
[i**B3] m2(D:^^i^V)\^7'-'inas»iVVy 

«ti-c*t), m2©;^'<i'h;i'f"-i'K:feW5StR3nft: 

tifctJ- :?'T L' figt?* S C i t SIS 

immA] mistitcmn)mic)mmtm2(Dm 

hMSt<!:H2©+7-b ^. hKS[i©ra©1iiE«:«J: 
3r^«f®/cJii)©m 

mm 6 ] St;? § n/ciJ-:/r u ©S!*sai^©«^w 

f :?'T U ©?^ J: 0 S l » C i ?r!^i; i -r 5 5 K 
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tc!i>(^:^^{>}-)Viy'y h'&^t^cti^tscti'^ 
mtfinmiiciimDZ-<i; h;i'7'-if^#T©/cs!) 
©*& 

!a©?SS^ST^-f7^-i?*«*&-r^;t»«:H>fiaf-x h 

*®ffl-rsx7^y7-i*^t?{:<h«:!RFmi-r^.i»^i 
fctatsox'^i' h;i'7'-3f^«TOfcto©*i4. 

?:S:»)©B$Ki L/T. m2©^ra?:f^i^®^rai L/T. 

mfiasJRsnfci^^tfffloi^HiL-cilSfjL. sjia* 

-5? i-MKffi^sn-sc i^fiFmi-raif^JS i ccia 
mm. 10] f-:^ hJi/T'-f *5, VSi^tir^ 

■9-7'TU-r©*©^Sii^-eiR«3tl. 

7^- i'*5fi^«®ti^5-«J; 1 § t»aiS®-9-:/lt^>®/c 

a>ccff i^y v snfcv f fx^T^ji^amm^m 

t-r •SiS3}c]19K:fBtg©;^'^i' \- jI't'- i'^m>tcib<D 

imm in ■y-:?'r b -(nfitti'^timwcomi^m 
^jj^etso. ^n-en©fjiBii^©*«c{3:^)£©sii[®it 

friaJt^«»»JSJ|g(cj:4;^'; ^ 
« (c J: «3 mm met i<m ttzm^m i o ccieis 

Biriaff;£©f»«a^B#ssii-c*»), 

50 *7^S«M{c, BiilB^«©ISA^^^SIJfC**ci{c|aH 



3 

ii©iifsmi©si^*^i©iii^K:ji^L> mDmiim 

2 ©^^^2 ©iti^tca^-r S XT- ::^<!: . 

t?C<!:*!|#itif-2.i«*:S 1 2(Ct^|g©X'^ii' l-Jl/'f- 
^f^^©/!:*©:^. 

Y, = a,Y.+ b (dY./dX.) +e 

Y,«Bflia-9-:?'ti^tc4BW5tii3';^'^i' h^br-i^t?* 
0. X,{JHflriE-9-:/Jt«'©4'©x-<d' l-;Ufi[gt?$)D> 
dY,/dX.«2iMi[-C*0, aSiyfbH. 
l:bb/a*i|«^W«li^:7 ^<&^■rJ:5^c>'^•^^-^»t?* 

(a) M!R3<ait1f:/TP'f©-€-ti-encc 
few 3 n/tiSilWMti-CctBlS-r ■5^iie«J-9- P'lf^ 
©:*cgf3?:«l^-r-SX7^-^7-i. (b) mm^:^^i' 

©■en-enKjPt or «jtaK^wti^^©;«ct 3 L -c 
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4 

SS? § n/i:8ina-9-::^T U ©^iT,-?txCC>i=t L ■ct^ria''^• 
^ - it?: n > - i'tt»-r 5 fcibtCmUmm^:^ ^ 

(cs^-r-sxT-i'^'t. (d) mmcmw%6^(ixk^ 
^iB^riE>'^■^^-^?Jt?^®^ur, BS*>6Biiiafti^© 
^^mm-r^7.7--j-p't. (e) Bffiatt->-7f--< > 
i'cctBiE-r &MiE3 n/ca;;*:^«t^^&ifci?)tcHiria(i^ 

^ffiffl-raxf^'^'i, (f) Xf--^^ (d) ©41© 
10 <fc'3/h3<«:.2>3;rmriaxf-2'^ (b) ~ (e) ^i^il 
i^rl^ffii-rSit^Jll 3{CiBtS©X'-;i' h^Uf'-.*^ 

*f©/cj&©:^. 

{4g©'en-en©cti{cx'<i' h^i/t-i^^nas-rsfc* 

17] K 7 yimtiK mvi^tifc-f-fr u 
20 ■^fig©-en^n©*©X'^i' h;naaKta{o:-r-saffj 

WH'^CX'^i? f- ;l/f-i'0 1 oCDi/'; -X^&l^-T'S/t 
«)K:«?rSfiT^*^*^^-C* ^ C i ^iRfSi r Sit 

1 cciaK©x'^i' f-ji'7'-^»^*f©fc&©:&!a. 
1 8 ] »*ijs«g*5> ■9-:7'r u 'f©*©x 

'^i' h;H4g{cS5f-rSX-xi' |-;K4a©fcJ?)©^>iftiKiE 
f±g*W1-S-';cX^i7 h;i,t--i5'©«?iE?!§€:^*. ^ 
©^^*f^iK:£f-,T 1 ocDf-i'X'^id' FJKigCC 

30 ss^ijsn. 

ifiai:;^^-^^' vs\^\^~^^7s^^ hJUfiigJctBiDSii- 

(a) mi©-!*-:?'TU'fi4g© 

'fag©*©frgaii'Xf-i'K:teWi->XX'<i? V)]^^ 
-fi^iRHb. (b) SSl©tBMSS(ao*-e©1f!£ 
§n/i:^J^?r«l5£L. (c) mfgafflKMI^. «fffB^3!l 

§n/cx-ti' h ;H4gRo-Ma«*«©w;eiitcj:«3 . 
i5?Kmi^*n-@L. B?iaKi^«;Jcifc(cs«>ifit»^'j>c 

40 St^jStRL. (d) friafflMMIK!. i[fB^S!l3n/cX 

'^i' ^•;^^4g. Ro'Bjgaafcifit^miXStKjco. m 
■rsili©fflBass4*tt#Lr. ju^jifeifit^m;?:^ 

( e ) 1 :^Wi/7 h itifcmm^i>'i5.\.^^W'm^M 

irxDWJi^^-imM^^L, (f) i5[-^© 
7;l/©*^*i!. tBB5Sii!©1^SSnTl,>2.^Jg?K:StL-r 

50 •^m.f^^mtl:^7-'^-:f {e) ^m.{^. (?) »]© 



5 

i^ntc^mi^-i^:^^{'i-mmv:.t&mb. (h) m 

ctiicj:io^-^m<Dm^m'bL. (j) mam 
iis?:ffjfi£-rs*rifie;^f -v^- (g) . (h) rj>- 

(k) MffitaMMI^> friBWJ^BH 
:x - i?itSL-rs c i *#mi-r €.ff*^ 1 7 ccEiX©;;^ 

T^^*lll©SaK«I^UTt,»SP^. lll©-^Xf 

-i'&g. 0i©^jaiWSi{RaftHi£;TS>xis?:Ji^ 

(k) ^jl^U. 

iria^ 1 ©->>cb-- i/fig. HiriBii i ©w^taMSf^s 

m 1 ©(4g>& 3 > f :» - $H-gt L , 
HuaT^S^-«:®2©aSfCi3I^Lri,iSffl. 112©- 
^>:f-i'fig> m2©#^WSMRi>'ffltE-rS'^fS* 
Itig-r €.fc«>{c4ift;®5T^(c J: ( a ) - 

(k) 

Bines 2 ©-;x eg. (jsas 2 ©wjaitsMSits 

S2©iil«:n>ta-*JH-^L. tJiem2©{i««JiB1f 
:/TU-Yfig{Cfi«?^t, 

*n-en©tsjt.-r-5f-:^T u -<(4g{cjpf Lrn 1 ©bjt 

$ffi4-e©tBJE&-r5ll2©<Bi©ffl©fiilli. Mi2»l 

©-^j>:f - i'{4giiriE02 ©-;xf - i?^m.tm<o 

*©i^s©^b«. ■en-en©)itR$nfci*-7'7u-f© 
*©-eh-en©siR$hfc-;xf-i'{c48its;^-^i' h 

JK4g?:fS-C3 > f * - $lt»-r 4©(C<^ffi-r s /c* 
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[it*JB2 0] MtC, 

Tg^fg^ ©f^oiaS{ctBit.t *ffis©a*? 5 n/cB^w 

{C. tBlt.-rSf^i^©-^Xf-i?fl[g?rSfeSffe^©-;^ 
;^'<i'h;bf'-5»4jR*ll/. miiBffe^©-ycf-i'<ag 
<!:mi©-;Xf-^{4g<!:©fffl©f^©M. ®Jr$©ffl 

-5/c»KMfBfSiK©M i ^n^-nciurhi © 

ti/£S5iE1f:/T U -Y©«©'en^ti©Sf)?3 nfcr^X f 

*{c. iiiE$n/£:)ijT*fcj:i5;^7^".'7- (k) *j^-r 

<!;«:^tf C t^fifffltil-^ft*]® 1 9 (CIB^©;^-^ 
h;l/f=-i?»*T©fc&©:^, 
Ci#^2 1 ] ^»j^$|g©*©3lc^T^T{c*j 

jf>^fcjf)©:^-r* o r > 

ift/Tu-r^rWo. ^n-en©H>iiat'-:/ru-r«f)iB«i 
x'^i' hji'f-i'© 1 -3©'>i;-x?:ii^L, Bifia^3fe 

a'JSilg«iifa-!f:/TU-f®tt3©x-^^ H;H4g{c>i?f-r 
S;^--!^ h;l/{4g©ft:«)©4^ft;«?iE^WL/. fu^^^tJ 

^^gwMcci^so^nfcx'^^ h;t'fl[g*wr.5-;>:>^ 

'^i' h;l'f-d'©SiEy-::^*^i?>. ^n-en©r;^f 
30 -i7{JtBH:£|g!itfr^^^*^^-©§f^JT$<!:K:«eor 
tBMKSiCciO f-i'X'^i/ h;l/{±g(ciSI»iJ34ift:SfB[ 
-m^m-i-i). x--;^; F;Ufilrg**esi)SfcJf)©7jffi{c*j 

HglBr;^x--ii' h;i'f-i'*HJIB;=^-^^ f-^bagcctBit. 
S-a-S/c»fC*^*s. (a) ml©•t^:/T^'ffil[g 
©4i©-^>c;^'^i'h;uf'-3'. 'Si.ifmt^KO'^^rv 
-/(4gso'^©ffi©si'?§nfc-!f :/7 u tta©4i©r. 

40 L.. (c) mmmm^.mam^iti^c:K^9V)\^ 
so'iriBSfeifii^sa'>:!e({c J: o . tBic-rss i ©fflw 

ialSl©tBMS»*^6)iS2.ij-5£?^*l|iiU, (e) 

1 h $n/ci(nBgfcfic>!g3S'X»%BfiS-rs 
50 rxowtji^^-nwm^m^b. u) ^-^mom 



t?> mamj: 1 miiyy b :^tiftWtJ:^^d:micJ: *) 
;^7-2'7' (e) ?r^igL, (g) m(01-o(DmV{$tl 

;^-<i'h;l'f-i74^S(JL. (h) iriBtBHB9IS!, if 

«riBtiH!£»{c J: 0 > Htc-ecDffe©^^IBl4 3 > f * - 
tmi/r. •f^JJDW^clK-Stm^rjSJlt. ( i ) #S 

-^ss^©Bi^ias^€riffi/J^u. ( j ) frian^cDBuia 
TS*-ciriB;^7^y7' (g) . (.h) Rtjf ( i ) mM 

(k) BijfBtBMMtiC. HulBW^fflMSij. 

n/ciri3i»-:/7 u ^ (Dipcy-^ti-^tiomsEmiRS titer. 

■r-2>/ci*«:iflSilrl*iaJf${cj:f5;^f--^:?- (d) ~ 
(k) ?r||3SL. 

iflBH 1 ©-iX f - i'&g> BulBll 1 ©W58!lt@MS»R 

cxm-r &^x^4r^ffl bxm^micj: o iriBiijf^© 

0 1 ©fig* 3 > tf a - ^?H-fi[ L . 
i^ia^^m^^?rm2©saccfil^L.r^,»■5.ffl. 112®- 

(k) ^H^b. 

friem 2 ©-^^ If - i' &g. iria^ 2 ©w^tsM^tSfeS 
cKtBiD-r ^miimm bxmmmmic^ 0 BtriaaJr*© 

»2©iii?r3>fa-$ft»L/> |fria»2©fi«i[iB1^ 

^Tu-ffiig{cte#L, 

■en-e'n©m-rsi^7'r u-ziigfcj* Lrm 1 ©®»t 

©-^^f-i'{4giiriail2©-;Xh'-i7fiirgi®ra© 

fig^i. ^ti^niriBfii^iKrtaagMtc^ifci^: 

^©iem®^J:t«:, -eti-e*i®jitRSnfc-9-:/rU'l'© 

cfi®-en^n®®t;?sn/c— ^Xf-^{cfcit^:^'^5' h 

jUig^f^f :3 > b- ^ - $ltSLT ■2.©{c«effl-r 5ft:e*«: 
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n^mz 1 (cieis©;^'^^' h;nas?:^jf>.5fcit>©:^ 
©f^©aafctgr£;r sffi:t©sit(? 3 ti/cB^ra 

CC> tBft.^*f^^©-;Xf-5'14a*^i?)et^®-;X 
X'^i'hJbf^-^f^rjRmL/^ «riBf^O-^^t•-^7^4g 
tmi®-'Xf-^7figi®fflOffeKOM. SJr*©li 

E*i^-r 5/cft«:Mias^©s 4-en'en©mi © 
t^®!aaccTgfeigT-*«ij#-rsF5. ^n^n®s»?$ 

^a/ci^iB■^^:/T b ®*©^n^*i©SJR$ n/c~^Xf 

«)«c. ^iE2ti/t®JT*{c<fcf);^f-5'7- (k) *^-r 
i^r^tf C i*!^il[i-r.2.W3iaB2 2 KiBlK©^-^ 

[m-m. 2 4 ] ^fumi- 1 . m^^^^tf^ hoy^fs^it 
%sffr-5*^tasi4wrs^3fea«£^g4*<ii/> m 

20 IB«IHJSKaiS©ftltBl*-:?'7U'^*WL, -eti-eti©*? 

i3t*-7'7 u -/«iria^as{ctews-en-eng:fes(4g 
k:*^;^-^^' h;u»«T©fcS>©J|g{ctet»r, 
m\<mmc. iHR3nfc1f:/7U-r{4gK:teWSK'; 
7 hS?i®/t«)©l|l©X'>;i; hJbr-if^riR^-rS^ 

mi©^7 42 2r Hr-3r*f#S/ce6K^n^tl©jg5?§ 

nfcBfriBif:/7u^fctei:f6ir'b->-rj§S3nfcmg 
m\<ommc. a»?$nyti^:/7u-f{4gK:fet:tsK'; 

30 7 h^Ji®fcs!)®^2®X'^ii' l-;l'7'-5»€:JR«-r.5* 

n/tifiBf - :/7 u -f k: m s^fagi (jibii 2 © 5- 
Buien 1 o)msiK.itibx&M<DmiR^titcmrMici£M(D 

- f i ^fijffl-r i *:a{f -r ?>Ct^^^tt^:^ 
'^i' H;l/^ff©/ta)©^, 
40 [iS3jt^2 5] fiJffl-rs^S^J. «S©-9-^7U-f(4 

g«:*Jtt4:t7-fe7 h^SrSfeSm©:*- 7 
f#.5.fc*{CMl©:i-7-fe^- f-f=-5r4<^-r^m ffi 
,«©+>-:/7U-^{4gK:fcWS:*-7H2 5r hi&^fe.5»2© 
hBgS[4f#5/<:»K||2©:t742i> hf'-f* 

^•rs^s. Ro=x'<i' i';ui/7 h^f^s/cAfcfria 

lll©+7H2 5. MaiK<i:B«ia»2 0:i-7-b-^ HM^i® 
^2 4CcfaSg©:^'^i' hJi/^W®fci?>©^|go 
50 7hSig*5, ^i©;^'^i' h;i'7'-»{c:fews K'J7 h 



9 

7 H'*iBi^iB*IHJS^ffi©/J^S C^-gp^J^fiSTS C i ^ 
!^<b-rsft*312 4 (ciBig©X'^i' h;l'^^®fc«)© 

mm, 

[il*]S2 9] Sf?Sn)t-9-:?'7U^©^*5. ttl^© 

t6m^m2 6(ciim<o:^^i^ hji'^w©/t»©sig, 

[if*:S3 0 ] MK. MiRSii/cWKfx ^•9■>7• 
;Kc*^WS7^x^;^'^i'^;^7^-^^*ltS^»•r•5#©. s 

* C i ?:«f8S[ i t 2 4 (ClBtS© ^ h 

[it^S3 1 ] sec, 'pts.<^i,\ rxDmmtitc^ 

d'X-tfJl', StR$tifc^F^{C7^x ^7'-i^©/c^?)©f• 
x^:^•^i^^;^7=-^I^JRaTS^S, MIBSlR^tifc 

Ra^BfriasiRS n/t^WB^©^gs*a■r/^•9 ^ - ^f^tt 
*&-r*fc«t){CB«iB-f;^ f-x'^^ ^;^•r-$^c, iriBS/7 

mm 3 2] < i 4> 1 -ocomm tifc^mi^ 

r;i'*ff'5fc*(c, ||i©B$Fa*5^fiJ©^ffl-c*0. II 

2(omfsmm<Dmmvs>t), mtR^titcmrMHi'pmo 

-if i-i^tCisE^Sti^C i?:!^®<fr5if3}t^2 4{C 
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§a«S©X'<i' h;l'^ff©/fc8!)©«g. 

:f7H'©tt3©ffStt^-ciR«sn. ^«;^'^d' hii'f' 
-3?*5^^ii;9-K:l:bLr /h§(,»gil:©-9-?'ti^©fc» 
{cff 6n> h Stifcv h >j i'X*7';u*SHijlBiiffl;^ 
h;U7^-ifAi60f^lt^©fc&{C®fR3nfc*7';l/ 

it*j«3 2{cfBlS©;^-^i' h;i'»*T©;teat>©Sls. 

•SCi*#i!ti-rSi»^^3 3K:ia«©X--«^ 
©fci&cSlgo 

[ mm 3 5 ] ^IK©-9- 'J h <t f) 

$ n , BifS©^*sa^«#siSfc* ») , 
«:*aH-r-5;^-^i' H jHia©«i->7f- ^ 

S. RO'BuS§fS©I^A5^^^-C* S C i (CjPf^ SliiE 
^/c«?>{CH5fBW->7r ^ >i'^ffiffl-rS#a?::i:fl 
•r*c<!:*!Kfii<frSif«:^3 4(ciBiS©::^-^5' 

20 ff©/c&©^g. 

m.mm-:fv^ifihWki i oo-> y -^©Kria^ft 

iJ-::^ii^©*n^nK:tew.5tfr4>oraiR3ti/cli««: 
©fi6©a»«irt3-y7'if^©*iB(jt!j?r*L. «riBJi«j» 

30 

fceftcc 1 -3©ia3gfctt:^'<i' l-;K4g{cHLr^n-en 

©HuiB-f-i'V';-x^:i/7 vt^^mt. 
|-»^i«^iE^ctew*^©»2©«^!£:ll3l■r*fc^<){c. -e 

^©B?iBm I ©*^*ll 1 ©itJitc, ll©iriB»2 ©« 

40 BuiBipF 1 (Dmrnmrn^ 2 ©it^fctew sssiL-s^si 
mi^^t^m\^, firiB3^ii*». «JiB?4i'7f ^ 

50 [if*]l3 8] l-;i'7^-i'*!;^scc«j:0^$ 
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n. 

Yt=a.Y.+ b (dY,/dX,) +e 
Y,BiiiBf-:?'tt^{cteWi,it^;^'<f h)V'f~i'-C$> 

i>^mt. (b) ^»MIS[X'^;5' b;l'f'-5?^:iR«L> 

^iRft-rs^Si. (c) atRsnfcfriei^-T'TU'f 
'Md/ HibT^-i/i^ririEHistcjs^-r^^ei. (d) 

■C. iiD^*^e.B«IBft^©fflll>&Bi^.TSS^^S<!:. 

( e ) mm^i^y f- > ifiimfc^-ri>miEs*itcW::k 
mmt?,tcibicmEmmm-r6^mt. (f) 

;^7^y7' (d) ©iti©it«3Jp?)©i/>3!pJ5:S(@M«>. mh 

y:/ (b) - (e) Ctifcj; OS^^tt^jl 

3B^S^S<t?:^t?c i4^t-r«.if*a3 5 {CiBiS 
©::»;'^d' f-;b»W©/c»©J|g. 

iritfdfyiC'pU < i 4) 1 o©^*Ti^*^t?SiS-9->:/;l/ 
•C*SC i*#a4-rSit^2 4K:f2*g©X'^:i' hiU 
»*f©/£«)©SIS. 

hJi'f-d/©! -5©i/";-X*||3S-rs/ci!)tC7feiffi 
{clBiS©x-c^ h;i'^#?©fc«{>©^g. 

[»*]S4 1 ] ^ms&&m.ifi. -y-^/r u -i'©*©x 

h;H4gcc5Pt-r.5;^^i' h;i/ttg©itJ?)©^ff5^ 
HiriB^7l£?fJS«g«il«C^lJ3n/cX-ci'h;l/ 

i,tctbic:^^i> f-JU^^fW/cfeO^gJliiMCc, ( a ) 

H 1 ©•y-:7'r u -f fi[g©*©-;Jc f- 5'{ctews-;xx 

^j, h;Ur-^f. ROfBUlBUl ©-l^T/TU-^ttgRtX-e 



C7) ^$13^ 1 0- 1 42 054 

12 

i.t?>r:X:^^^V)\^'f-'itt^n.mh^WLt. (b) 
B?fBM3tifcy^--ii'h;K4g. mifcor 
m. 3 nfcSfi!J©tBBII^Riy''?-©i|$^ § n/cMJiK J: 

•)i*iiw;^i&*tf-»^i.^is<!:> (c) 

«:gfejSt,»^;>:i(?rSfR-r*^Si. (d) S4>i5 
l»Slgt;^:S!tRO'^ 1 ©fflMSI»*>6)!X*f(-SfgtM*ll 

m.immt^h^sn.^mcm.t^^ o tus-r -s.^ i ©ta 

HS^^rH-gt-r-S^ISi. (e) iyvVt1nfc.mJSL 
10 O'tSlE-r StBB8S^*^ 6fiXSl-:5©M)S'2.m - Sii>Ff * 

^•rsfc&fc. 1 h $nfcHfliBgfcfic»^is 

nn^mt. (f> is-)eiK*t©»^*itBM5£ssj© 
irsB!|$!£$nft:««j^©/c«)«:ji^$ns$T> iimmu 

1A:W'>7 hSnfcSJ&SSflStimcJ:^*® (e) 4 

iSiiL3iffl-r2.^i2i. (g) mmmtitcifzfr 

u-f-C. •aifi!i*JiEK<J:0^*6li/ii<2;ifii;f-i';^'<i' h 
;i/{4g(CBfl[iH^'>:x'^i' hii/f'-f ©— h;i'f 
-d'*i6»(|-r€.#®i. (h) (=«jnwfe^-^fg:>pf* 

T?$)o-C, iffBtBMMISt. B0fB^&ft:f-i';^-^d' F;Ki 
SSiyfH>riBSI^©^<i^ti©tBH®K*Mffl-r€.#l8 

i. (i) ^¥m^-vAi-^micm&-r?>mii<^mmtj: 

K$L-C. ^-SS»f©«^«:ffi/ht-S^I2i, (j) 

Hfria«^©*t©*©/cfci-3©tBMSsgi*i^-o-c. BiriB 

30 ifc^mig). (h) Rz^ (i) 'kmiKL^m-Ti 
*e<t> (k) -en-e^iosjRsnfcifr^Tu-Y©* 

©^*i^n©H9fiBStl? 3 nfc-^X b-- i'KijWS X-^ 47 

h;i/fiig*tt«f'5f^®-c*-9r, itriBtiMMij. friB 

-;^b--4'*!!. ag?r^-rX'^;^h;Uttg?:W 
:^^i> hJl'^W©fc2&©Kg*5M{c. 
40 tBIc.-r SlPr 1 ©-'^Cb- f (ag?:^s!)S^ 1 (D-^:^^ 
{> h ;l/f=- i?, 15 1 ©W^WSI^OftB[D:-r^;Xf^€: 

KfiBmi©-'Xb-i;{4g. mumKDmnm^mLVf 

©H 1 ©fig^r 3 > b :x - 5?H-fiL-r i . 

tBit;r.2.»i©-y:b-d'fig, m2©w^ffiM^S 
i>'tlit;T«.iX^*S*S^2©-;X:^'<f hJUr-i?^ 

50 mi^-r^tcibi^mKDmLiic^mM^^mi^-r^mic^ 
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immmicj:*)^^ (a) - (k) il^bmt^^ 

t>'fflit;-r i mt&'AWiic i 0 )BSf 2 ©ffi* 3 > f * 
-fitSfr-?.*s-cs>or> g[ii2offittBuia-t>-7'ru 

n(Dm 1 ©sJr^ffl i ^©tKt- 5» 2 ©fit i ©ffl© fii 

friE^ 1 <D-A b--i';^'<i'h;l'fil[g<t iriB® 2 © 

n©iiiM <!: (is^ i omti^ > f ^ - if n-g-r 

*©*si^©MJ:b%> •en-e'n©jit>?$n/cH!riBif/ru 
-Y©*©^n-en©jS!)?$ nfc-^x f- i' «:*f-r -s;^-^ 
i'h;i'ttg^S-cH-SLT.5©ccffiffl-r-S/cJf){c> ^:/T 

mkt-r6m^m4 i (ciatg©^-:^? h;b^#f©fc&© 
[i««3«4 3) Mk:. 

g<t^l©-^t:--i'fi[gi©K©Mi^©M. fJIBf^igS 
©^©aRC/^n^*i©^tk*=i>fa"ifif-ffL/> s 

*^-r S fc«>«:^<t©«EH-Som®«f n > f a 

* nf'n©siR § txfcmmv-:f7 u -t ©*©-en-en© 

Sra(C4:;f*®l)T${Cj;0^S (a) ~ (k) 
©;X'<i7 h;u»ff©fc*©^go 

ji>*fc«{)©sig-c*or. 

:/ru-^%Wi/. ^neti©if5B1^:?'ru-/«i[iB^ta 
Sfctei^ 5-en^*i«ij© 1 o©{4g(c*o > HiriBT^s 

Ttt. MB^t!f|g^ROt!riB«ltli»cJ:»3. BolBt^^T 

^'^©t^l©x•^^; h;ufii[g{cffit£;-r-s>^isijfl«)-^;cx^i/ 

h;l'f-l'*>efiX-5 1 o©->>;-X?:j^-rS/<:«?>K* 

is^sfiu. fria^jfea'jssgtJ, miw:rr\^^<D^ 
(DT.^^ h;HafB{c*t-rsx^i' hJKig©/fcjsi)©^:fift; 

^n-en©irferiXf-i'«, tsM^iSRo'B9 

iET^^^©^^Jf$K:fi£orffiMMtS!K i 0 b - 
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(a) rnKD-t^TV-msOtpomB-'-Xii-fiC 
is»i>-:X:^^i'h}i7'-it. S5faill©l^:/TU-Yfi 

gRo'-e©ffe©jiw § n:fc'9-?'T u {ag©*©";^ f 

- i' {c*j W S - ^ H ;U7'- $ *iRII-r i . 

( b ) friBtsraMit. mi s nfcfr^;^ h ;i/{4 
g. mi>-o-cmM$iT,ftMmomm^m.^'t(D^m$ 

10 miism.m-^xmi^mi>m>mi.^>?:msit^^Sit. 
( d ) Ah&^mi:x^Rufmis^ 1 ©taMSis*^6 
^sis-5£St?st**s^s;rc*(CBirfBtaMM3S. ifsn 
St$> a^iJ3nfcBiriB;^'^i'hJUfigRofg*>ifii,»irfa 
mwi'^xwiic^t) . tiis-r i ©tsMSf^sifmi-** 
Si , ( e ) H snfcttiia^xi^Riyftaic^-r ^fiiB 

itK.. 1 /cwi^^ h $iitcmi>'M>mimm^xwL^mi& 

i> (f) WL-mmM<DWt^*i, mmmwi<o<wmsti 

s^Si. (g) ;3ij©jiJ»?3n/fcirfa-y-:^ru'f-c. 

Ma^^ft;KiE«:J:f)3}c«{)P./i:4i?*f-4'y^-^i' hJKJg 
(cifian^xx^^ h;i'7^-5f©i:^;^-^i' hjuf-^;?: 
liSiJ-rS^Si. (h) #ftlW^cS5!-^SyFf^l^-r 

■C, HulBtBMMiK. H5iBii;^;^-?i'h;K4g&CXifiBll 

^©■eti-enctBMse^^rfijffl-rs^mi. (i) # 
bsj-c* s ^xSKtBic-r s ^=t*nw^c>pf ©*©i--< t:©ffl 

^im'hiri,^mt. (j) iri^^©*©BiriBM© 
fpcofctc. 1 ':><Dm^mm-:>r. miS^iE-^:f7 U ■< 

Rc;-e-n-en©3i5? $ nfcmria-!^- b {ctett swjsj 

tBMSISROftBlt^S^O ©gSi!^XI0[4:J&fi)E1-^*-C. 

¥13 (g). (h) (i) 4i^taffl-r&¥iR 
i. (k) *n-en©sj?$n/cfriBi?-:^ru -/©>*» 
©^n^n©Ufar;^Xf-i'{cfeWS;^'<i' i-jnag* 

40 iai^5SWBM^Riyffii^^©H>iiB«0©SfB!'XI»%W 
ffl-rs¥Si%Mffi-rsci?:«®[i-rs;^'^i' 

g*5}?«)-2>fc&©$^g. 

X-^i; h;U{4g?:*»^fc«)©£?g*5H{C, 

m^thn 1 i'ffig. m 1 mmm'm.R 

m%o^m.mgxmK.i.*)^^ (a) ~ (k) ^r^iguii 
50 ffl-rs^si. 
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muf^ 1 f-i'fig. mmm i ommmm^R 

taiD-r an 2 cD-yc b- - i'fig. m 2 ©w^fflws^s 
iy^fflrt;-r.2>fiS'j>Ji*s&-sn2©-;x;^'^i' h;i/r- 

T^KKirie^^iB^JT^KJ:!?^© (a) ~ (k) 

«5ie» 2 fig, mriB» 2 ©wjsifflBajsssjR 

$©m2©<ii?:3>b-*-i'tfgi-r-S#S-C*-Dt:, ^ 

-e-n^n©tBi£-r -2. if :/T u figfc te w a ^n^n© 
mam 1 ©fflff^ffl 4 ■eotact.-r s Msaii 2 ©jB i ©m 

2©-y<:t-i'x^i' h;H4gi©ia©{iiaM. Rif^e 

^(omm<Dmt^, -en^n©siR§n/cHiiiEif:/Tu 
-Y©*©^n^n©jijR§n/i:Hir^2©f-i?{ctet:f 20 
siitB;^-4i' h ;i'fl[g?:Sjritgt-rs©(c«iffl-r5/c* 

il$ffi[<!:-rsff^:B4 4{ciBK©:^^^ h;i'{4g***.5 
S»f$©liiE4JI^-r i/cJ{){Ci[fBfSI^©->>c f - i'tt 

Bi s 1 ©-;x f - <4g i ©fa©f^i^©M. atFiBf^i^ 
©^©aRiyf^tx^n©^i:b*3>f :x"5fH-»L, s 30 

^ n^-nojitR § n/cfifB-y t/r u ©*©^ti^n© 

a-i^n-m-rs/cfttc. f§i^©iag-cTS^^^^4i»ij#^ 
i>mc^^mm^ici:n^WL (a) ~ (k) ^m&^m 
mt^^Wit^^tsc timktti>m^mA 5 ccib^ 

(07.^^ h;l'fig45}?®-5A:«>©J®. 
[it3}aS4 7] i^mk^t. miZ9>fm=i-tfth<o^^ 40 
3t«rSftf **^ffi|gi?:^t?»^feffl;£SlgO*©-9->^ 

fc&on > b- J. - Djtg)feiaigj«mctel,»r > 
iiiBt^fflg»a?^©t^lB-!f:/TU"f5:WL. ^^i^ti© 
ir^tf :/T u-f «ifia«^aiS«:fcW€>-en-e*a«iJ© 1 o 
OfigKAO. ^jfetJ^ggfJHfc. mi©^raccjllR 
1^ :/7 VA^mcmZ) K 'J 7 f-«tS©y!:»©ll 1 ©X 
'^:i'Ml'f=-5'4iR»-ra#l3. m2©B$W{Cjiil?-9-7' 
T u u -^agfctew.^ K y 7 h^«i©fca6©l^ 

2©X'<i7 f-;l,7"-$;£iR^-ra^©, i)iaif>:/;l'* 50 
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mmi-t^:^'^^ h;Ut»$ll?:3>b-*--3»|+g1-Sfc 
ist)©;^-^^ h;U'f-i»*S^-rS3>b-*-^jH-»^e 
iriBiBtijif*«. 3>b'a-5fH-gt^®{Cj:0 

^ 0 nfffi-c* Sfc«)CcfjiBiattm©*cc^j!fe$iTC 

- $ a - K «> ^n^n©s!R $ tifcmttLV^r/r u-r {c 
*jw.5i5«>-3r}is§nfc^fi[g%^*. BfriB'?"c3i^7 

^fimR s nfcSifiBit :rT u fete w Sir -3 ■cj§;£ 3 
n/c^figtMian 1 ©x'^i' hMT'-ii^tmt^^ 
mt. 112 ©5^:7^2 v hf'-i'^fi-sfc&K: -en-en© 
if >; > ymriB-y-r/r w -f ^c;tew■SHufB^^4g*H^riB 
ii2ox'^i/ h;i/7=-ifijt$$-r-s^Si. Buianio 

{igtcJt-r'S;^'^^' h^Vi^:? h4?f4A:*ccHuiBlll © 

ffl^ <!: C i ?r«FSS; i -r a 3 > b- ^ - $KBS{ 
OoJtBi&fBtmfto 

S©f-:?'7U-fttgKStt-Sm2©:t742-^ hf^-$>& 
m<^-iVM'j |-M»iMfBll2©:*-:7H2ir hMIS:.!:© 

pa©^4Mffl-r s^g^^tf c i -r ^11^4 

7 (CiaiS© n > b-:x " ^fKKO Bltl^iiaiSjift. 

[iS5j?]B4 9] si«5nfci^ia*5mi©i^Him2© 
^Pfli©H{c4.i3. lli©:*-:?-^ 

e>n S C -r 4 7 K:^©a > b- :i - 

if^HS{9pJ«Bi&:iBKi®i*. 

7'ciy7A3- K*5M«:. tJIBf-x hX'^i' H;U7'-if 
%lSffi»*jMSSIg^ff={c^li{k1-S/'c*«:x-ii; h 

i'^ h=&ffiffl-rS^g:*dtfC<!:?:#Siti-rait*54 
7 {ClBiS© =1 > b- :x - i»RK f) Bl«l^j:iBttiift. 

7=- if 3 - K*!gK: . ::j>^i: < i 1 ocWitR 3 n/c^«f 

m9(Dm&imt^^'7 ^ - •5i^«>«:friB7'x 
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[ 5 2 ] ^*iiJ«gS*s. 1*- u ©*©;^ 

msi,^mmm.tim^. (a) mi©t^:^TU'fffi 

1 ©f-:/TU'f&gSD''^©ffi®jlJR3n)ri:1J-:/r u-f 
{i[g©*©r;Xf-i'K)i*-r'5— ^Xx-Ni' hJi/f'-i't 

(b) \mm^. mufitcx-<^v)\'^m.. mhr, 
xws. § tifcsfij©ffiM;tiKRo^*-e-©^$s 3 fifcw$,^ 
jd'omm'^m^^n-irh^^t. u) m^mm 20 

( d ) irf as «>aosi!{;>:sisouiam 1 ©tSMSi^ 

mi>h)S.i> 1 o©gfj:SfK-S|g[S=f4ieiTS:fcai){c, 

<f) ^-%m<o$k^ifi. nm.m.oym'&ttifc^ 30 
is©yt*(cji^$ns*-c, #]i)ne«j^ci:r£w*/7 
fcm^^mkm.j;:nm%i^WL (e) frigiiLSffl-r 

S^ISi. (h) W»j:»-)£IS*if*ll^-rSfc«> 

latBMMiK. S5ia5j«;t-*'x-ii'h;i/fii[asc;frie« 
^©■en-en©tBHSigt**ijffl-r**s<!:. (i) 
^fabs;xftcctaio-rsffln«%Biiia*f©4'©-r'^r 4o 

7;u®»^*^eT'<r©ia!-SIS>Pt?:BiJ|^Lr. t^a^ 

-s?a*t*ffi/h-rs#s<t, (j) wmm^m^i^ 

<Otp<Dtctc. 1 o©taHSiaASJ^-3 r. ItriBfiJiEtJ-t/T \y 

Ro'**x^n©jii»?§nfciriai*-?'T u kj* lt* 

i. (k) ^nf'4i©SJR$ti/i:BirfeiJ-:?'ru-r©* 
0-e*x^ti©®JR5*ifczi;;:b--i'{cfeWSX'-ti' h;U 

^. iffBW^ffiSfKSiyfBgia«^©H5fa^o©^lS[^>: so 
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[ if S?« 5 3 ] i . fiia»f$(^i^*> 6 o^tij 

JK4g*^** it«!){cffjffl-r S /cad© 3 > f a - 

mimtiimmmo^tiiv- :?'t u-y *w 0 . ^ti-en© 
©{4gcc*f), BuiaT^^iJ. Huia^t^J^^somria 

«J|(C. fflM5£S^iB<rfa^gN^^-OW^»T$4KSfo 

H>iia^a5S^g/&SM{C. (a) lllO-9->^TL"f& 
g©ttJ©->Xb--d;{C*jl^S-^^X'-ci7 h^l'f'-ifi. 
mfsBH 1 ©1^ :7'T U (4gR!yt-?-©ffe©jill? 3 

* s {c if lai Btsjif*© * (cn^te 3 n-c t i S -r - if 3 
- KRi>*:/n A n - K . miBf^- n - K 
«. B!iai»-:7'TU'f©*©x^i' h;U&aK>(rt-Sx^ 
h;H4g©it«)©frfB^»J^^g©^^ft;feJiE> Ufa 

BJiaT-n y 7 A 3 - K*i. ( b ) mamm^, m\ 

3 n/cBifiax h ;K4g, W -3 r «e£ 3 nfetiTiBft 

ii^mi>^mt, (c) B?iE^flf)^xigtK:g4>ia(,^ 
^i^c^^sfK-rs^si. (d) Mciii&i^m.m. 

Rl^mim 1 ©tBM^I^*'6SEaiS-^l^?r^^-rS 

tc^icmmmmmm. m^mm, ^iJ3ti/cX'<i' h 

;KiigRt>'S4>ifit,»SSS['^IS{c i 9 . tBfS-rsH 1 ©tB 

m^mnnti>^&it. (e) '>7i-3*i;fcffis;^ 
i^c>m-r4«riafflM;£fBi*>^fiSs 1 -5©M!&:sif[- 

iK^wtf-Sfrs^Si. (f) i^-siK[)i*oii^*5. 
«iiatHMSi^©«FS 3 titc^m(cKi bxmm^ti^t 

Mia^S (e) *iigiiL/®fflT«.^S<!:. (e) MO 

/t4>ft;h--i'X-^i' h^HigCCiriar^^X'^ii' hJb'?- 
^<0-JK7.-<i;V)\'\^~i;^mn^^WLt. (h) 

ft;f-;:?X'^i' hJi'fiigSi:miB«^®^*i^n®tBM 
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ti^mt. (j) miim^(D'P(OmBn<Dfp(Dtcfc 

(g) , (h) Ruf (.i) ^mmLmmr?>^mi. 

(k) ^tl^tl(DmVi$*lfcmii-f^7U^(D'P<0^ 10 

S i *$tf C i ^ri^m i f 3 > f J. - ^BX 0 pltfeife 

[000 1 ] *^Bj(j. :3-»JS^g. i Ofc)l:f CCDJ: 

«)CCipJfflT-S3>f:x-aj^f)W^&iaffiiSftiKM 20 
[0002] 

[000 3] 1 o© f -/©^J^KSfl^S^gtt. M-?©X 
^i' h;u^*J^3i-rsfcii?)(c-9->7";uaA4W-rsg|i» 30 
ie^:?'5X-7 ( I CP) =&jpjm-r€.. I cpiffl^#r 

otctbomamt. 7,^i> h;i'iSi©^a'J3ns<imo» 

©gS7feifii^*wrsi?fmiS^iis (cdd) ^©if* 

m.i)imiini>. C©J:^^ctS^=&iHji-Cl»S3X^tS 
h D^-^'«^R@i|ffi=S4 8 2 0 0 4 8#M 
iffl* (B a r n a r d) tC^7i^Stiri>S. 40 

[0004] fses**j«L.< uitmc^mmms. 
®«©^». so^ (0!iA«figR{>*E;^i^tti*^i=)©) -e- 
n-en©^»j5t»s©*© k 7 1- ^ciigsi aj: 

or^fc. *li|#«F»5303 1 6 5#MifflS {Gan 
z^© "GanzltitWfflS" ) «. aiK©^«'JSSI 

gcc^ji© 1 -:>©it5t3n/ci^7*o7 

{Cj;S^«JSS|g©^ii<t4ia7Kl/t:i,>S. iJ«J« 
«g«cj: O^sns;^-^^ l-;u««IIBR{cttWffi©li 
RC;f7'D7^;U4W0. ^ss^btt. :»«J;£^g:?"n 7 50 
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^;l'©i^3©^KiK:>tt-rSffliE*tf ^. C©J: -J^r^illfk 

h ;l'©ll^©jSEfittii:'-r L fe!t:<Sr«ifel^ G a n z 1^ 

Ganz!f$S?giffl»». igg«iE©fc 
*©*gHiJ WiSS f - i' © y - ;^ ©ffim *> L r I » 
L/*iLC©J:^!feffiffl«, H2yy>h{t$n;^ctJ-:/TU 

[0005] •»^>7';^©4l©fiS^©^¥aso'JI«. x-^ 
i' h;U'f-i?^©^iEt-f;i/©jSffl«:J:i)3>fa-5f 

yfcMfi^snr^ l^>:7";l'©'t>©^iJ?:S*W{cgS(ifi<j 
CC#*, ■?-n6©a?r^*'2>fc&©-!>->:/;l'7'-i?CCji 
ffiSnS. lo©r7"O-5^«*a#«F05 3O89 8 

2-^?gffl» (I va 1 d i^) (cga^sn. C©r:?'D 
-5^(3:. •»f>:/;l'X'<j'hJl'7'-i'©^M^?:vn;i' 

•S. ctilix^i' h;b-r-$*5, *7';i/©i|J©^Mia!© 
x-^i- h;l/7="$©jR;git^lJS!73ttii^;StcJ; 

g^^fbs-ttscisciojaisfifjfcx'^d' h;i/%iii^(c: 
-5f^|Rm■r■5©^c^lJffl3ns^c^>*»*:>b6■r. xy^y 

[0006] 

■5^3t»ISSg©*-C©X'<i' h;l'7'-a':9'*f©Ai»© 

*©C©J:^«:7^Siyf^l8*Stt-raciK*5. S'J 
Ktf'ibSEf'- if «r^«Ejitf WfClBS-r •5/c*©C © 
«> :^»J5^g©/c4!>©1?JiE*^Jt'?rBfiS-r*fcii?)© 
Sam. :^ii'J5£SSg©«-C©X'J V VM^iimt 
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[0007] wfjiisma. ^7£tJSS^©*-ciJ->:/ 

[000 8 ] mm<OR:/^<Df^(Omit'ptj:< 1 4)35^ 
[000 9] miOl^^i' H;l''f-i?*i»l©B#fflCC, 

{ciR«§n. •en'en©sti?$tiy!:if:?'7u-/K:3«-s 
K'j? h«iii©fc*ccjDm3ti. ^n^n©sj/?3nfc 

» 1 ©^ffltc*f L-Cffi;g©jiJ?$ ti/c^KttEO-y-:/ 
WfiJ{C{jmi©:i-7-fef h-f-i/it. ffiS©f- 

ffiS©-9-:/TU-/&S«C*t-r-5:t7-fe»^ l-=Sr5£«) 
Sir 2 ©:*■ 7 12 5» h MS5[*f#-5 s^cAfci^ 3 n. 0 1 © 
:i-7i2 7 ^H^5[i^2©:^•7•fe^; hMfgti©^©^*. 

[0 0 1 0 ] i':? KSf-X hf'-iftcjifflStiSCiiai 

;Kciiffl3n-5. c©Be;j©fc«t>ccs*v i- 

mV{$htcmiiiC7-:^ h-!J->7';u©fc«)«ciR«3n. -e 
©^, ;<'^d7h;Ui^7 ^?r^t/r»*:x■^i'^•;^■?*- 
jgffl^n. C*iKJ:0. ®ll?$tx/c^W^©iSS=S:a-r 
[00 1 1 ] vh'j^x-=e7';l/©fc&©T-*-f:/T- 50 
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7£tJ^i^g©*t?©x >; ■)> h^SKj: oiRSisn. cn 
n. ^©^s 1 -3®iii#{cto/c9^;^-rf :7'©raKgt 

[0012] h'<}y mm$tifci^m<oit^ 
T-riScife*^. -e^TenotJ-t/Tu-Y©*©;^-^ 

©«jEv-;^^^tf. -eti-en©-;J:b--i'«, (^gt 
S-?-©/l$fcJitT.5) #-©iaK€?^cctJf->r 1 o©f 

/tAW:— hJl'T'-ifjJs. l|l©-9-:/TU-r{4B 
©*©— 'Xf-4'(C)ipfL/-rjR*$ti. hibf' 
1 ©t^-yr ^•Yfi^gscK*©fl6©s^'i^t^/c■!^• 
:?'T^-^{^a©4'©-'Xt•-:?(cs^L■r^R«3ns. n 

w> ^-s?^©«^*i. iii©i>-7'Tu-r(csturti 
MS^©!^5e3nfti^i^fc*H/-ci^snsst?. -wtt 

s. 

[ 0 0 1 3 ] »ij© 1 o©3is?$nfcif :/r n'ri:;j:f 
-i-*!. <&ft*^iEK:J:0^«)6tlfc<&^i;t•-d'^4g{c^ 

msti^. mmmm. ^i^f-^^fig. so^jt^^©-?- 
^■c*-s;>cs;K:fflit;-rs^W'S>Pt©4'©-r'<r©tB 
ll^©«©>i=t©4i©/c/clo©tBMSIfc*s«-3r. «^)) 

iisns. tBMMft. w^tasa!£igt5a^«'^©«©«f) 
n:fc-t>-:/r ©*©-en-eti©ji!R$nfc-;Xf- 

K«. fflMMfi©*-C^3n.5T*^-?©TO**3S 
[00 14] 



mmo^^mmmmi o mi) Kiowfflsn-s. 

^3feaiSS?gBil^;ltf, Perkin-Elmer O 

ptima 3000 ^©K^-xti-f t>mm>&^r 

S/c:Sf)CC7fe?g©fc«?)©-9->7";l/aA«:W-rSiS*iilS^7' lo 

mm(om<i^pfm5 3 o 3 1 e 5-^g3iffl» (g a n 
z ) cciftH^§nt:i».5 4? 7"3?©^^ffi©/c»tc3i 

i»5^«g^b©/c»K^fflsnrc>s, 

coo 1 5] C©J;^ife^g©t|J©felfflSl 6 (02) 

wam« (*a[SA©fra©*s^s=»4 820048 

^^^m ( B a r n a r d ) {cpg^StiTteO^H^fflS 

(cgifflK:j;o-e©^»*jaxjisti-ct,>2,fe©4-r-2.) , 

©•b i'^>h{bSn/£:1f:/TH'2 6iLr©2^^ffi 
CC»»Eg3tir(,^.5j^3tttii^2 2*Wr*> iEL^ 

WK:«?^me^ilg (CCD) C©^©;^'^^ 
i' I- l^;ili-St©tS^XB 1 -3©^^ i 1 -:3© 

[0016] *Wffl»Ryfi»*©«H©*r© "■»>-:/T 
U-^" <!:©^«« »iJ©li^{C«igJ^XB2^^«fflK: 30 

2 4 5©-9-:^7U-r©*K:7 0 0 0<DmMi^tsCbi> 

io. •en-en©i^-:7'7u-f«i o-4ooiii^$:^ 

(0 1 ) L r f# ?) n:*:»»J^i^g© 3 > f :x - 

3 2©>*';H2i'i^3>3 occtSiWsnrt^s:^'*^' I- 40 

;l/!Rfit2 3©fc«?)©;^-^i' hJU7^-$2 1 (03) B, 
l^^^Jm?S(0Wil^'^:f7U-^ 2 6<D<p(D^tl^tl(D 

jfjs^gBi/EioB. mm&micKiLxmm (xam 

^) K*fL-T«E$n-2>o i/*>Lf^-r Sixteen 
B« *^»JSSg*sjS;S7'-$*ii:>^i-t±-rtcffl^tff 
il^rtlP.-r ■5*:%BJ©t,. < oA:>©BS-CBi£J.S-C:&l^ 

«FStB'>!5c<<!:^>fll*HCliSS[©ai^{Ct>/cO (X«;^-< 50 
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i'h;nag©*Rgti^©*-c) i^yht^ (%mmtj:i 

fflriafflSti-ci^s^^fcB "jj^g" sex "jS^" i© 
ffl^isjiis©)S^«j^x«^«*;si»-r^)&^. X 

t3*^aSTb-1'X«-9-:/TU-fffi*© (miiEi^:!- 

h©*©) x-^^ h;H4g^3iB*-r.5>. ) . 
[0 0 17] n>f3.-3J3 2BjimB, mBD i g 
ital Equipment Corporaton 
model 5 1 0 0©j{-*tcJ:0*^3-»J^B 
©•fKiffl3i**T.r(,iSD i g i t a 1 Equipme 
nt Corporatonmodel 5 10 0^© 
fie*©*)©-C4)S. 3>fa-5ft3:. (MS©J:^{crt 
)asC>'/X«i1-#W©) SW^cT:^ai'/f=^yj/jU^ 
^**-r^.**ill»^g (CPU) 32. T'^Xi'* 
^tfCi4>*S>*';42i't'3>3 0> :<--^U-5»A:^ 
©/c«i>©=^-!j<- K . RO^*- 4? - 3 eSOVXBf^^ 

- $ XB*^»JS^g© ^-*KJ:iJ*3>fcr*-i» 
©tf«:ji»«iHiiStvCi>.5 "C* * " ??KJ:S*>©© 
J:^k:fie*©fe©f*.S<, *%?^'s©7"Di'^5>i^© 

[0018] eEai©HS3*i-Ct,JS7cmfiS©Jg^^ 

tj*^i'j^g©*cc^A5*is. mti I c p 1 4 

©* KSA 3 It, 5 2 0 . C ©SW©4>©7ES©ia^©a 
fi<J«i*ftl'5!£>g«&t,>. /c/dOcn6©aWffl3trS®WAd: 

H^-T^B. t^aiS-r©CiT.6©11^^7C^©X'N;i7h;U 

rStt©^«*^?&3nsi^{cjafR3nrt,ie. lo© 
(cJ:f3:S'-y-^ hK:$tiri<>*-y-:''TU-f{4g3 8©S 
■en©jiti? $nft:-yyTu-/©x ffi^so^ y issfc ^ o 
[0019] mwjimmitm^ii. 1 70x^9 y)im 

tji\.>Cth$>i> (mB2 4 5©^{CJtL/rS|5^^^ 
©tfifC 17). 

[0 02 0] vu-^^-v (05) ccfcc^r-en^n 
©•y- :/T u «ffii*XB^{4g 4 8 ^sii^r 6 nri* 

fp. <KS7'-i'tiji:Stt. -en«n©ii»{c*fLre:t 

-5c<!:«>$>«j, XB-en^n©ru'r{cfcw*x'^:5' F 
;^!RF^a©/c©Ji5fiS<agso•Je^^©a»WI5:x»•?•©ffi© 

B^tC. t?II^«1^yTH'©*5^l?©-y-:/TU-<©ffi:t 
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[0 02 1 } x-<^h)iT-^4o ms) omiom 

^'-^mWi^SU. ii;^«i«#^Cj;0^{4g4 8<!:4 6 
CCffiffL. Olir«1t-:/r U>f(Z)xyJ^:S3 8 (04) 

lCD:t7-fe':^ h7'-^:*^6 5 4'C5Rj6?)n (n>t'ju-- 
hO = ax+by + cCDJi5'C^O. 

SV\tfj:lK 30 
[0 0 2 2} F'J:7 h^*iaf 'S/cifttcX'^^ ^;^f'- 
vy hmm4 4<Dfcitf^4 2vmm^ri6o mzomm 

^mm(^c4ev^&m4 8tim$ti. cti^j:^m 

2(Dty^yhy'-'^32f)mf^ti^. ty-^yhT"- 
^6 2t3:. m2<D^miCi6lii>i£M<D7U^i<LStfCi6lrf 
^ty^v h^SZi?)^02(D:t7-fe-> hMSi6 4^5 4 

hffiiS:5 2. 6 4CDr5CD^6 6{i. C(0^ y ^-^^)V(0 

[0 02 3 } StRSn/cB#H ( *'B$ffl3'* ) SOiCfeli 
xmk^tifc («M9^J^c*^CD) hlf>:/;U7 8CD 
X-^i' hJbT'-^?? 6&i. V)\^i^y VI 2(D1 4 

t?£D®ffl^C J: 0 1 4 3 ^(DftltCD-^--;^ 

B$F5ccswb2n^ci;«ii8 2ccfci>rpiee-c5>'2»o ffl so 
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^cmionmtm2(Dnf^t<Dmxt£^m\tf^hf^^\ 
c(o^^y'<DWmit\t. S5<D55i^^c<fc*:)^3nrc^-^ 

[0024} ®iift«W*tJtCJi. '^lyy')i'v'-^{cmm 

*B^i^^05 3 0 8 9 S2^mnm ( I V a 1 d i 
"I va I d ) J. C. I v a 1 d i RO^ 

T. W. Barnard* ''^^WT'-^mmM^m 

(Spectrochimica Acta, 48 
B. 1265 (1993)) <b CC^Snr 1>-2>*;S-C 

[ 0 0 2 5 } V h U ^xtJ. rX >:7*;KD4iicm-:> 

WCC Djfier^-2>o I V a 1 d i ^?Ptiv h y ^7X(D*CC 
i-:>tLhcD«Mifi:^^i?)'^Ci4>IS^Lru-5o L3&^b 

^) o 

[002 6} 

imi] 
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[M] = 



ai6 



^16 ^Ifl 



[R] = {[M]^[M])-«[[M]^ 



[U] 



u, 



[002 7 ] [C] - [R] [U] 

[0 02 8] C(0«^m 6©H^&Cj*L.rT*-5. 
[0 02 9 ] (mf^(D I V a 1 d i ^f^<DmmW)m 

c<Dimittm:^mt^. mnmmm^ms 2 1 s 5 5 s 

(de Loos-Vollebreg t^) 
Sl>'E. H. van VeenSO'M. T. C. de 
Loos~VI lebregt* "il^^:/^XvJ[i 

;l'^f«;>y" (AnaK Chem. 63 144 1 

(19 9 1) Ruco^Xifcmm^tixi^^^Miomi^ 
[0 03 0 ] ^mno i -ocDW^fj^csg^ ) -ctiv 

ccvhy ^x^7';i'8 4«. -y-^t/jbr-^? ecDiRI! 
m^hfc^^^^ h)i'^y h 7 2iit. hSti/cvh y 

i^X^f^vTb 8 8 ^mmt^fti^itC^ h y ^:K^'fML 8 
(B$ffl3TJR^$n/c) ■tf>:7';l/7=^-^7 6tc9 Ori® 
tJLf'-^ X7'U-r9 2ra^Sn5c J5iJOlo0rXh 
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[003 1] y d^;^^f^;l/?:g3|fe-r^^. 

h7&^ (^ftfl^JtcgcO T-:^^^:/iR^7x--Xra(cKy 

[0 03 2] vhy^x^f'Jl' (07. Wm>7.'rv':f 

WRiWa (^pfli) 4 3^7&>oiijR$ti/cru^tt 

g3 8 (04) -CKy? h^lg4 4CD/ces?)0[)*7J*3>^^^7 
30 h;bf=- ^96*42 riRHf -5 C i CC<i: 0 ggif^-T 
SJSBX-:^' (lya^O) ^{4g4 8i4 6-e 

-5^7-fe-> h?:^i?)-5fi«!:f7-fe^>. hKIStl 0 0^54 
rff^ (a>t*^-ifH-ff-r€)) /ci?){c«fl33n^o ? 
ffl;^-i^h;l/r-^l 0 2;!^^;j^C^t?*^B$Pfl (W^S) 
80-Ciit^?^#f^l O6(Dg£^0rSffi(cS*UTJR«Sn 

-^10 67&^. Ky^ heW4 4K:*fLrf$^BSrfl (B$ 
40 H2) 6 0(C4 2-CiR«Stz. f^iKO;^^^ H^'I't'-^ 
ti^{ig4 8i4 6t?tbl5$ti. cncci:0f^CD^7 
•fe-y hX-^4^b;l/f^--^rl 0 8;«>iff6n. fgigSC[)^7-b: 

'^^^'imO^V'^iv hM^l 1 0^5 4rff-2» (n> 

5?tt-SL$n/c:^MSS[6 8tJW^7 0^C<±:0?iJffl§n. C 
n(cj;0ttiP5B#ffl8 OJc*f-r'£>*PiX'^i? h;U^7 V l 

50 ii^f'-^? 1 0 2cc7 4raffl5n. cncccfcDr-;^>r 



[0 03 3] fraiCDj; ^KttJffl;^-;^' h;Ui/:7 ^^:f^:@ 
■CAS. *7';l'?:f#S/t«)©a7©Hirj^x^y:/*iS5 

-7h ';i';^*-f;i'©iSffl©B#ra(C5^7 l-i&ig*cc3>f 

yit#©©H{cl2«©:^7^ :/sy'**S»c©f<i^* 

[0034] <fco^<©^«?ai 1 6*i#?i-r.5*^ (a 
mt^nWiEti) i©«|igi 1 5K:ic;i;-cc©:?"ai' 

•ji';^! 2 0©fcaii)©T-*-^:/ 20 

■f - 4f©iRft(is$?7-r s r c n ho)^m^Ru^m(o 
7 hifi^fii-CB. ^oi^'jy hifi:rvby l-mmJ^Kifh 
n^n©^M5©^cc h'vyh mm^mmr s c i 

S3R< . K :? hS*i*i'©gg©«arffiffl-ri.3!)^K: 

i'h;i'7'-5r©-en-en©ii^*. F'J7 h^ittcjro 30 

[0 03 5] 7-*^7'7'-i»4*«)€)Ci*s*J7f •5 

( 1 1 8 ) i> Si|i?itii«*5. fii|i?liffl©fc«)«: K >;7 h 

•;i'X*7';n 2 0K:)PfL/-ctfton.5. f^SWfcftfg© 
X-^i? h;!/*^? h©^f|:*. f§© K 7 biiJ5£«i*>e>© 

©if^©ga©«ig<!:L-ctf5ci*5^fie-t?**. -r 40 
[0 03 6 ] x-ii* V)\^f-nt. m^tm<m%:^^ 

i'h;l.Jt^©*©li^fig(cJ:f)tgM3ti-5o jgfi (X 
S<!:WipJ-C*0> C©i®^, *»»JS»g©St^ffl:t» 

htt, ^liSR©c|3©7»-$U'<;l/tctel:fS/h3<,» 50 
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mttbxmii>. vh'Ji';^*f=Jl/©fc*©»*;^'>< 
h;bf^-i»©T-*'C:?'f'-5»JRSI©fc»{c«J:OiS5 

i^!isg*5Mi;ns. c©@e«]©/cJ?)tc*»7Ea'j^sig(3: 
ffwtctj "igsgg" ^fis-cffiasn. c© -mnc m 

h^^-r -5 c t ccfc 15 ^pjter * 0 . Alt;^ ^ M 

7 (01 ) BXf f f>d/t-:Jf 1 2 5{cJ:i3*i:a*](c 

^waif'p©foii?)ccji-r s.-i-rapi-cpiset s, 

[0 0 3 7 ] asCCtel^r-tf^iUAS, X'^i' h;U7=-5' 
^f-X HT-Sfc&Cca^SiaS (9 0) i> ^y=)\,<om 
h;l/7'-i»8 4«B?ai©<t^{ci/7 ^§t^S 

(8 8). il^»{CteWS*'f;l'7=-^?rffi 
ffl-rs/cl:fCJ:l^ fJ^b^, J:0Wc>53'»?li©*7';i' 

©«© 4 §a«©f= - i» ^ > h ^ffiffl-r s/cw-c i 

-fiji. ^n^n©efi£J9^©*©a» (mu4^) 
©if^ccssfLrff en (■et-C'V^ i-^i^T^jwr-^-f 

:/iU«©ra©Ky:7r--f >i^{c*fLrMESti) . 

[0038] ^l^©7-x h-tf->:/;i/K:>Ptur^ 
a'jsns-r-<r©pjegft^*?^. T^siyws^r^tfx 

'^i' h;I/©;^#(,if'-^''^-x-r«:to^T-*'f :/*siU 

ftsns. cneox-^i' F;^^^:!;}!^ ^:*csNj:b*s 

K.isi.^XitmiRStitc'm) ©«Se«J«C«6 0~7 O© 

f^^wmM § nr . v f- y 5- x ^•f jp© a:»© 7 - 

T-*-f7'©/£»{CffifflSnS<i^ 
©^M)1^>:?'.'I'«. 3iS«frai©J;^(c(,K-5*>©^}> 

[0 0 3 9 ] ■r'<r©7-5!7-^-'^■;^x•^}' FMi. 10 
©fta©»*sjssig^«:a^iiMsns c i AiBiie 
■c*s. m^immici*)^ty-ky F^WTssms* 

?:W-rs-r'^r©x-^i' fjw, ^^r©x^i' F;i'*i 

ia?&^{c*st>r-k>5rU>ysti/ciieg'^.i/:7 F§n 
•5). f^i#fltJCc*>■oWfU{CBT-;<7'f-'^*;^X'^:d' FjI/RO' 
tajt-f •S:i-7-fei>F*. JRSiS ti/ct»JS-C*Siffi^ S C i 
*5Dlfie-C4)f). S*3n-5ffit©i^7r^>i/^. 5> 
$ AK:^&6n/c:t 7 -Ir h ( F y 7 F ) ^^{cff -tt 
-SciJtiSplffir**. ■e©fi£i*{cio-ca:S^iH5«. ^> 
*f-r'^^S^©'4J©t'^r©fiS:»-©X-^i' F;U?:iESI 
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MMfMftJkit. CtiibOZ^i^ hil'Ot^XHK K';7 

[0 04 0] mmicj:^:iry-bv hiiiEtt. miiA. 

S a V i t z k ySO'M. J . E . G o 1 a y*©l63S: 
'mftsnfcg/jNi^:;^ni'-y f ic J; S t'- 

tfifRlfMrnW (Analytical Ch em i 10 
stry 36. 1 6 2 7 (1 9 6 4^7^) (Cit^S 

nri^SSa V i t z ky-Go 1 a yfi^JI5l?©fie*S 

m^mbrmmstii. 4-^>*->'^•;^;^y•:» f-^a 
is> 4#@s©.#.iS:iijRffig(ctaiigTS> -r^ito^m 

©*> m5©*> if9©.'^.lf-C*S. ^©ffi©«M:t7" 

[oo4i]Rmmtj:}=')yvmm 20 
2 (08) ^ii3i-rsfcj6©7i£^Taf«T-cAf]. cn 

fctt{?J^tf > -e*i^ti©-9-:7'r U-r 2 6 ©4J©fi[g«:^ 
tyx-ii' h;i'^«{ct)fcD*iIiJWf£:raPiJ£gi,^r(4g-r 
SffilS©T*SSX«-'XSfJb-- ^W1-■5lil>'^•$ 

"(57 ^4^;^x5fn>" ©}^©$S«ji©^ 
«, HirKE©*ai^?Fm5 3 0 3 1 6 5-^?a*i» ( "Ga 
nz^¥fmmm" ) icm^tSti. C©J;'5«cx3ia>(c 
M<St*S|J»«*?«iffi»©«{c?|ffl(cJ:i)Wji*n-Ci<» 30 

[0042] ^tl^tKOHA^^- i'W. f - i'jffigtcM 

orssstjpise^cssiK©;xi!(*«-r-6. mwj:mmmm 
{cEssnrt,* 3 - f- > y $ nfci^ * 

h-C&S. ( "7^^;^ (F i n e s s e) " 

HS©jifiiiefiKfflPi?:gti-rEg3nT:i»s^ 40 

nri^S. ) :t7'5^*JbhL'-> (01) 
|^7^5'j<>hX«-9->:/Jl'3ftlLO:^^Xvl 8*^6© 

I- 1 7K:JR^$n.5o lEL-t-^n 'J>-J^a >©fc*(C 
:i-7'5^*;uhu->©i^©xifn>l 2 4 (09) «2 

■o©i:t,j(ciaffi©u>xi 2 6. 1 2 8©r^cc-9->hv 

•^/g^KBBSSn. C©8f5*AWU>XiX'J y hi©ffl 
(CieSStimS. IP^l©|!ab>Xl 2 6-C5feiSBW 
(CStl. ifn>*Sjil/, ®2©CiU>Xl 2 8{CJ: 50 
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[0 04 3] ffi7^*Xi$n>«5fe«%. «l^©SSf 
-i'*JfMT•ST^^^^•^'->©*K:eiT^. c©^ 
e^^^^-'^>^-7••;>^'©ai&^e*>x«G. Ghos 

h. M. EndoRiyfT. I wa s a k i ^ "(,»< o*i 

-<^ti[SDffe^>f5(" (Lichtwave Techn 
ologylS, 1 2. 1 3 38. 1 9 94^8^) tcgi 

w § tit I* s -fe ;i/ V >f ■Y-~im^m<Dmm-fi^f)^f^m 

[0044] mk^x^mirAommmism^icj: *)xii 
n>js t {cjcSD-c-en-en©^- i'a?sS[{c*fLr^ij$ 
US, 

[0045] 

o = m/ [2n (cT. T) t ] SCl a 

tctd bUmm^ a RCJffigT'NCSJT^n ©fl!c???r^ 

T-S^tISggO*©x ^ D >©ffl*>^fi2[fil^ D 
«-r6. J:«)-«5W{CB. Ganz#i^?§iffflSOffli5=S: 

[0046] *:^3fci'JSSg«{S7 ^ ^X ( "L F " ) 
j:ifa> 1 2 4-C{^S(lSn (01O©A). cn{Cj;»3 
ry^x-^i' hJi^-^ 1 2 9*iff p>n> ry^X'^i' 
f'-i' 1 2 9ttilS ( ) f-i'^^l/. ISf-i; 

iii^K:*tLr5i823n-Cl»^.f-^{ieTl2 2 (0 

8) 'kmtio ^1 a«, m«mD©^*it>->:?";n 

3 3. *^«JS^a©4'©RFg!Iii^>7'. X«09© 
:i-7"9^;«7;l'hU->©*©fi7^*Xx3?D>l 24i£ 

g^^srtSKftg^i* (0iiA«G a n zmmmmm 

^SnruSNd : YAGife^) *if>©e{»©^l^®iR 

i^ii©gE»i©f-i' 1 3 0 mt'tom'C^mm) ©its 

CceESD©iSfi{C>CtL/-C|Hl1iK:4 2 Ktet,i-cf#?.n/c-?- 

^fiiffl^-tt-rftjffisn. cn{cj:o. ^flsn/t (tt 

m «*WT'2.MBgCc^a)6<ifc-^^x-:i'h;l/t- 
1 3 0?:a-r-;^:X'^^ h;l'7=-f 1 3 2*i^fiS3tl 

— ;^:f-i'«-9-7'rU'f{4g©5^©i-^r&o. 

C©-!f7'TU-<{4g«*ffii;-7'TU'f i^cO. 

B. S!>?3*lfc1^7'TU'rfig©^^©lo-C*So f 

-i^y-xttfeJIEv-xiLrttf^U. c*ik:J:I3. ^ 

s. 

[0 04 7] C©,'iS$fB7'DS^-i; + «ffj^©Gan 
zi^rFB^*iS{cl2tS©7'nS/-yf fca(ECC4>D> C© 
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SCX-eOiemOf-i'iggO/cJftcDG a n z 
[0 04 8] c©J:^^^f^#©7'ai^-i^•1'■c«> vftlK 

:t--''«7 -;':/-rSS*)jfil»«Sf-i' 1 3 4*J 1 3 6K:te 

(-^r^D^n. «iif-i'(cj*Lriaf-:?©iE?i?t{4 
©gi5»©iE?i?c(±ta*^«).2>fcis?)fcji-^sn.5. i^f- 

i'«:*jW-5iSS®COg|5^<iHlf«;XISm, + 5mK* 20 

m»;^^©a}»^t?abS. SCI aB^XSCtc^jrS. 
[0049] 

m, + 5m=2n (a,, T) t (To ^Ib 
x$n>©;f§ t.©i 3 8K:tei>r©WWJiS«. W 
t6^cJf§^i^©fca6©!^§n/cffigMJi-cffton. c 

ntc j: i3fi]«S/IW§ t = t . - 5 t *5f#6n.S. fiJW 

5m4) (fiiffl*:>e.) 1 4 0{ctei,i-c«is$n. cntcj: 

^■C^cti ) BSm, ' ©tte*i 1 4 2 CCteli-C Wte t ft 30 

TO-rsji^ t*i^i a3!)>6i 4 efcteiirimsn, 

Cnccj:»)»l©t. m,*f*il 4 7(Cte<,»-Cff6n 
S. fft^©Sf^m,*il 4 4(Ctel,»r -e© 
BS, ffllS-rs/f $ t t = t . + 5 t S-ctf-St 

$^•5. mtfC©jg?!g[R!yfffleT=2 6-C-C«?ifil!§ 
*afce^'J*{Cte(,jrCTo=3. 770739lxl0'c 
m-\ n=l. 50023 ■C*»3. 5t=10/imCC 
r. /i3«t.= 500/im4«lffi$n, fi]»i5m = 
0. 2 5<imiJi^$n€). yc(,i-Ct=4 9 0{C>FtLr 40 
m,' =5 5 43. 5 8 6-r$)0, ^©^. m.i>o5.<0) 
IgKm, = 554 4Rc;>[tlS-rS t = 490. 0366 
VS>i. ;>:©SI^«m, = 5 54 5-C*i3> ^<Dm. t 
= 490. 1 2 5ii+St3<l'S. cn«t = 5 1 0 (ifi 

8S[ • JI3 (m,, t ) *t©«^*s 1 4 7 K*Jl»t:f#e>4l 

m 1 ©ii 3 i 0 ^Ji 3 tcs* L r >; X h -r s c 
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jimf- :7'^HB 6 ©^*5 3 a i'as-Cft (, > J: "3 so 

[005 0] 1 0©B) /cfetcglj©-!?- 

:/7 u-f c©S«©/c*KSiJ©«?g{cte 

ft<. fiifS^aKiEXli^H-. X««aijL«et*o©»» 
> :7-;i.^©^|»*:. e. S ^ c i *i nJt^r * a o 
1 -:)©Bi^(CJt^-SiESItt-C3t^-C$)S. 1 4 8 iCfeC^ 

rSf-f-r^Tu-Yccfet^rn^^ (iS) f^-i? i 5 0*ifie 

tl, iE32ffliK*lC©7='-5r{C3i^$«i. ^t-i' 1 5 2 
*i> ;2rf*«iEX«t>-7'rU'f©*©gt*Df-i'*ie,^6 

tii,'<mai^tzmm<D'pvi 5 4{ctet,»-rstRsn 

5. *f©^l©y:^ h35>^©JI3fflt*i. IS^ffitt© 

n^s^^f^^m<DmmjHv^mti>\:^-'rAmmi 1 5 e 

{ctet>-rthg-rs/t*{c<jOT34iS. io©iii^©;FiS 
^14*. ®^©^KitmT*ci*spJtE-c*-S. Sl©/h 

fc) ^3nfc/i3*>e.«iS3tifcmffl©fii»-r-<r. 
3Sfi©7-x ^^ii®■r*fc^*«:ttSI©c©/J^3^,»^^i 

[005 1 ] ^mX^m, ©**iELl^O-C. #SI?g{ 1 
6 0-i)i^Wloni^.^t^tcii>icmii>hti (15 8) . 
SlOOt. m,*fl 6 4tt:*Ci}i{Cjlfel>3n2.. ( "^NSS 
««S©9T5tK;H;fC*4. ) c©7'a-fex 
B. fcfc* 1 o©5Et*5a-3t:t6n4 (16 8) Sri 

ti(l6 6)> cntcifJW^xtJfa^JlS t (X«ffl 
B8^S[2 t) j&JJBlESnS, ^m^ti>-Anm,1fi-^X\,^X 

[0052] X $n >©E©if 3 mmti'Bi.umM 
(DMWimimhtix\.^i,im.'i.xhi)^^jt>^c.tifi 
wtg-cftu. ?J{:9©/^•^y-^^i&iiE^i-c*ai^5S^^^, 

-i^3>^©*^^©*©^StStt*s, iS^?rfil?!PK: 
^^•5. J¥3***eK©nIfi6ftfifci>ft^ll«r&t?C© 
J: ^ ft^?tStt©ffliE:^«. mmo^ 5 etifc 
JlSti. n. = n, + ant?4>-SJ:^{CiaiSSniS»K 

sfij«:its3n/cs»f^n.ic©s»f*{c*t-rs 

©i^ftffllE^SIifS-r&^ffiB. gtftI©»iS?:Wr.2.^i> 

iK©(tJrae«j#«ti«:ffi«L. cn6©s©'&ti-en{c>(*L 
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m*iizmicmatim^m9m (uv) m.m<D(pvmh 
[005 3] x^ayisttcitrnmiv. miR$tifc^m 

JW^gJfeL/cJi^. K 'J :? HiiEeD:/n lo 
«. jHR$nfc^«f^©K»i®rSmc»L-c^6ti, ^ 20 

!K«:j;sfiiHKj:f)«t7]©B#raK:s^7 h sn-ciisn 
■5, ^Xi^-C'?^;^ f-lf>:/;u{cteWSf'-^?fRli©B!rRy 
^«c^iSi:/ll2©;^'<i' Hji/f'-if^Jii/n^K';^ 

*i. •7hyi7;^*7';U'f-5'*5tlit.Lr7";^h'f-f© 

*^?g©c©iite^!iit?Hx-^i7 hJiT'-^amm 30 
m. (Xim&) imm^ctifimtLi.K com 

[0054] mSME 

^t?>::^^{' hny-$©«©^fb{cm-CSgfe. 

T^j:fo%x3?a>{c*jWSaS©«l^fiIi. S 
SOM^iLr©SJf^6CfcW.5ffliE©ga<!:-C*S. 40 

©2-5©U'-c;U, 0^;S.{fi:i,UC2*ClSint:t^.S2-:>©U 
x5rQ>ttilS;S:^{bLi''<:*>4>Ln^Jt,»©-C^©* 

liCttl^xm^it 1 /i mt(rt©jfifi^«:fig-r S C 4 
ti,) . ffljT$©^{bKfe*^*t>6-rSSr$«l^ffln„ 
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[0 05 5] )K:)^©aS{ctei^r7"n'>-i^-f tt. 
©«gK»-r St"- 4^ *JR^^^ C i CC J; n . KS* Ur 

©«©*fl6-r.6'f-5?) 3!H«til3n. ChK:J:0> cn 
e> 2 -3©fiS«:*jW 6:^^i>V)V f - i'©SMSafj^i^ 
©<ieTK:st-rs?gfi©B{iifci Lr^isf^3)j»3n 

e.tf c©u Vi^ajagiStcx if n> b- i^KSturi/ 

ctxcc J: 0 )lJr^©^{t*iiaS^{tK*fr *iaft©Mii 
*©f-i'©;^'-ii' f-Jbi^? F-©^{Ccfc9^3*i> 

tct)fcsJiJr^o^{fc=&«-r43ie©«^ mmm ifi 
•^x\.^xm-^<mm^u "ij--*^-^?" uifimminK.m 

tbxmis.$tii,o 

[005 6] 1 -o<D':faiy~v>^ mil) *s. mmm 
■s. =F^s^*ii 1 ©ss 1 7 0 fciffijt-rsra. m i o 

OAtmi 0©BiSl aicmbXWmStitczfaiy- 
V>-^ifi. 0T^©^ft;)lSf$n.KJ:f) 1 7 2{Cfel>r|| 

if. mi<Dm}mmmt^RmBitt^-A^ia,i}im 
mtii. i«T$n,©^i®fii«, i(i©wsafflM!t 

<5 1 7 4fc*Jt,irrJ>fjL-i?it»Sti, Ill©«l«. 
•5. T*^^?:IP^2©SSCCiffi^L■Ct>-5ffl. 010© 

ji'f'-if. ^2©*^tBM;ei^2 t,som-rS'Xig( 

m,*JS|^$ti*. S»T^n,®m2©iaAJ. II2©W^ 

7 4(ctet,>T:3>fa-i?it»34a. »2©«i«i^->^T 
L'-<(4g{c«»-r5, ^E-n-en©»i©)ffljf^fiii^© 
m-rsm2©ffii©ra©<ii^n.-n,A5. -en^n© 
^m:.thV-:f7u^iiLmi,cMbxi 7 ecctet^rn^f 
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^&gM©^ltR= (n,-n,) / (a.- cr,) *il 
©Baii[i Lt:-e®IS^©^i:k«^J!i5 1 8 2 (Ctel^rlSfi^ 
[0 05 7 ) t'-**^. 3i;i;«*SD©safctet,^rxi' 
->>:>^'^f hJi'T'-fAn 8 4cctet,>riR«3n. cn 

(CctO. fflj£-rsfSiK©-'^f-i'fi[g0,*sS«?)6n 10 
■5o ^©— ;Xf-i^figilll©-;^b--i'ti[g<!:© 
P^©^©Mct3 - a,B, aS^<b©Rai U-C 1 8 6 

©tait-rSMtbl 8 2<!:©a^n= (a,-cr.) *R*i 
1 8 8{ctei,»r3>fa-5?it»Sh. c:n{cJ:0®*T 

$©?iE*J|^$n> ^{b<!:4^t(5HJf$i©i^H-n,+ 
^nJ!>n 9 0Ktot,Jr3>fa-iftfg[3n> 19 1K: 

ttiE3n/cS5f^»^>ci»-c. SI o©ASi:; 20 

-f©ipoSJR3tii^cr;i(:f-i'K:)isti/-cx-^i' h;Hig 

MB8»©ctirwfflsn5. 
[ 0 0 5 8 ] ffi:*: 

iIj$©J: ^v::k*)v Vi^m\t. iim©li3§7J1SK: J; 0 
nJtB-C* S *©»<i: O ^ < ©*§tt*S-r -5 jTc&tcSS C 
i,^. Ct^«^c*>^■r<. vH;i7;^*7';l/©fca)©T- 

fx h-f-4?Bjim«7;U©iiilgt|:3-®*-CiR*3*a 30 

m2©x^i' h;l/f-ii'J>. +Bj£-r-St*-:/r U'f©*^ 
§f!£©-f>f-^fJl' (m{3Bi«) -CJR^Sn. iHiX 
^i; h;Uf-$(a:. ^S©lt^<i:f3/jN3l*1iit©-!^yii 

[ 0 0 5 9 ] X h^Kfet^TAWX y 2- H 2 3 
(01) «Xf -^b->i^^r-5f 1 2 5CCj;f)SI:^|fiIKX 
y ^/hSn. cnccJrfjX'^Ji' lo©if:/ 40 
Tb-^cc** io©ii^©3li;?$tifc-y-:/Ji^ (/hSc^ 
^ffl) t^^«lIl^S4o©-y7'Jt»?<:ttKi*>§n 

•s. *Me«j«ci?fl©^s?:> t^^tfHssnytxyi- h 

i' f 48ii*>-r c il?©-e©ffe©^l8:{c J: ») ^fiSf c i 
ifia}^v$>?>, •r'^r©il«{ctew* i-:>©7Ji'©x^ 
i' H;nJii:«;{a-e-*af'n©-y-:/iii»t?jR^sti, cne 

©X'^i' h;H3:f^Ccffl^-tten. cn{Cj;0^K©1^::^ 

^e=)«:KPi^ai,^-ciEg§n5*\ x«eEto©Y>f so 
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4? ©iijfWN^CCiEJtCCSlR 3 ti ^ C i ?:i£>S<!: -r 
^. T^jrto^X';-^ ^■^^•?■^t^©S©Efi6^c^A^Sg 

h y i/x*f ;W©4J©tt*E3n/tS*X'Sj' h 

[006 0] ->;^^J^i)m?Hi 1 -3©iiilS4lttl*5-r/c«) 
(CESSK X y 5- h © 4 0 X f - T-^Wr-Stl^, 1 -:>© 
ia3^S9 3-5©-y-:/ti^{c5stL-C. 1 OXf -^-/HPit? 
jR«$nfcf-5?ttiESItC0. 2 5H^U«:fin-CI,^ 

^. L*>Li^xf ABji^«. iia^aofii^cc^j&.s 
mmo^mmmxitA o + e ( £ ti/h$(,»^^«i) ^© 

xy 2; h^S=&Wrs, c©JS^. 1 0X7-2.:7-/cW^^ 
St$*ir(,i-S4o©^Saio©ii^?rSt)-r. fi^o-C 

4o©x y i» h^snfcx'^i' h^i'©ia^-&©^w 

x-^i;h;U's©iEiNc©^«rapg©#^-tt-c*S, M 

iKj&e.«ffl^©X'^i' -e-n-en©®^©*© 
0. 10, 2 0SCK3 o©g:fii«:i^$n5*>^-c* 

fcx-if h;l/©J;DiEti«^31Aiff6nS. 

[006 1 ] ^m-zfm^^mw.hmitm^^'mrr 
*5#^SE-r5. 1 -^oymmm mi 2) vitz v v v « 

tf^^X»-r>5'-''^•J^■C. tf!l^«a^{CtbL/T 1/2 

Xtti/sxai/i o©•y•ylt^x«-/>^'-''^'^^■c 
xy 2. hxf '^7'©IS«, ^ft;P|j>&{iB*> 

li, mii4o©-y-7'it9'©it>©ji^©xy h^-c 
w 1 ass f) 4 0 © -r> if — '^';^©sl^©^i^*i^t^-*s, 

^S©^®ia(Kil>TjiJ?$n.5.Ci*il 9 2{Ctet,>r 

/^3t^^ill«ii©«^A5ji^R^t^. c©«^«-ii©^ 

f^©1f:?'ie^>S:5£».5, C©afB[ ( "-T^t^-:; i'X»" 
tmmtxi) tt. ^:^^^^;iliH•«:i:bl/r»(cfet,^■cJ:^3 
m«3 0-5 Ofc^i. 
c©^oitio^n-eti©i' y -x©^Stt2 o©iiso 
tf^Xf - 7-SK:tft>n-5. rfet>% 1 iSSSO 3 0 , 
3 2, 34- • • 5 o*t?©•y•::''-f>^?->'^';^•T?^)-^., 
[0062] ■y-:?'T V -r©x y 1 9 3 H, -y > 
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^mm^ti^. m^mo. szmmt^) ^tti^ 

(^ICDf^-^ 1 94 a. ||2CD7^--5?1 9 4b. ^3 
Ot-'^ 1 9 4 c^) i^f&t^. ^tl^tlOtV'Jt 

)\^<DT-^iyV-y<U. 1 9 6^*50iriocoili^/ct:t 

-So c<Dv'-^{t:x^^v. 1 9 stcfci^-c-eotBis-rs 

h$ti/c7^-^ ;E>>6«St5ti^. M 

1 9 9E3:|lfS&C^<Z)iS-CMet^n^7l>\ 
JI^(D^(D«^-S«]^:^ffi ( r s s ) om-C^^ibti 

7*© (mEi*^J4 0<D) ^^CDSS:tCia':J< i r s s 
tiii^-r-So 0iI^«IIKCD»3;^^iEffi^4O©ii^. r s 

[0 0 6 3 } it'^x T V y^ommo mmtimmo) 

y (x + ^x) = y (x) +^x 

C0067] 

Yt = aY, + b (dYo+dXo) +e 

^ h )\^mLm L rsf/c{ciiJS3 n/c (gg f u ^ 30 

ff^lB^ "0" ) f-^tt»gP(D*(DJ:t3«l^ (xy 

llJ^a?)4)o (Y,^C>^-r-5>7"-^t3:. 7-:^'Y>'^' 40 
^ ©iii^P^ PitC ^ 7"- ^ r ^ e C <!: * 

[0 06 8 ] it^^ttih^W:i<^m^timtyy h so 
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T^m^Cj^^m^hfl. ffe*(DS«<&f*>f>f^':r^^Xj:0 

:^^cW>f^';^i7X|^ce:J:0*J?>6ti^o cn6©2o 
(OmtSX^^CWMomi^mi^. 2 0 2CCfcl^r^r s 

0 3tC:|dl^rili^^0CDX'j h^2CDj|^0i6:^^ 

[0064] tfes^;3ij(Diim. Mi^mt<km^^m.^ 

[0 06 5] mi&<oyjm\tnm)xhh^chti^i)-^t>h 
fflu. ^o^. '{k&iyyv^i^hxmmmc^:hmiLti^ 

[0 06 6] 



s;2 b 



(dy/dx) S2a 

L-CSi2 bCC<:fc(3=J>t':x-^ft-g$n. 1^:/rU'rt^ 

fcm-(ommtifc'^':fr\^>(xx<. m-^xm^"-^ 

l^y htiX'j Vlt'^T.T v^(i^±^^\cwmh. iC 
i^T^^b/aiCCDA^^iCDra^Cfi^^n^i. StR^ 

(rms) {J:. iELliffi:A:^^6C0fli^^^-ro CO 
[0 06 9 ] cn^fflJE-r^/cJ^^C. ffliE^tl/cJfc^^Oi 

iissn. m«i. oo*^e>fnt:i^'So x-^^h;b 

2b:;5>5SO*®ffi$n. j^-^lth / ^nT.'r vy'<D 
^gg^S (bisectional search)^ 

V hccMLrtii^>^ru-^tDxyj^S<DMmr*^o 

[0 07 0] f?*ffiCC« (013) (Xt^ 7 7*) 

2 0 4(DlSilfl*}:^tai. ffljc£;-r'5iSilWi^ft2 0 6<b 
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1 1 4i L'r2 1 2CCtet>riR«Sn> ^Mt5:2 1 0;!/5 

2 1 2«C*j(,»-Clt»Sn-5. 3ia©;^'^i' hJl'f'-iJY 

-€'nK:jtLr2 1 4cc48i>riR«3ni>. mwo=f-^ 
saf»MI^«s:2bK:2 1 6{c*}t,»ra^sti, cnfc 

;*:gf32 0 4SJ, 1 -:>©fl^{c 2 1 8 CCfcUrji^S 

"ft^" BWWfCB. ^(ni^^ZM.m--n^ ( r m 
s) 2 2 0*i3>f:j.-ifH^$n/£:iti^r*So C© 
rms*S?SfCOr> 1iiE3nfc5yc2 2 2 (1.00 

*^?)©^a!i) *<2 2 4{c*ji,>-cfiis§n.5). c4aa> ij- 
:^ii^>©i)K«S$ tift::^:* S 2 0 4 {c<fc 2 2 5 tcte 

2 0 6*5ff i/-ir>;^». KH2 2 SKtSOr 20 

rms*i> Bq4>->-CStR3nft:U'^;l'©tf{CA'5i;'CS 

WMtsti^ c©^j^.f*kf^©^i?wMtt2 0 8*i, mk 
©iiiH-2 3 2 i 2 3 0 (cfet^r^toens, 
[0 07 1 ] %w!^mt%(mwm.^^m\u,xmm:. 

*©fiffl©tfK:AS^{tRO^M*sS^(cBeHgt?$) 
fi£-:>r*^?§tt. aSff©»*©©HX«^ti^©« 
ii{c J:-,T©*$ilK3nS «>©iT2.. 
[|SIH©ISWj:SiW] 

[01] {cftffl S n e:9-7feite^©ftlDS0-C A 30 
[02] !S3t14iigS©1^-?'TU'f?r^L/0l©S|g©** 



* ©*IH1S©HU®?:^^W0-C*S„ 

[03] S2©«^aigl©lS^{cJ:«)Sfi3nfc7ex^;l' 



[134] 02©«^ttlg{ctel:fS(ig)l^;&Wr.5i«K3 

n/tnf u ©gi5^«^«r/T^-rBS«0-c* 

[05] *l6?gccJ;2,01©«il{ctB{S-rS3>b-^- 

^fffgt© 1 r>(i:,^j(im>7.T- V -fRu^mov u-^ 

[06] 05©7n-5^ + -Kcfflic;-r^;^'^d' h;i'!tt 

[07] *^?g«:«fcS0i©gg{ctii£;-rS3>fa- 
i'ffgt©^']® 1 -5©||J6JfJ^®x-f ■T'SO'^g©^ n 

[08 ] *^Bg© 1 ^©.igeifcjpjfflsn* K y 7 vmm 

*(=>©Jiaijfl«jX'^i' h;l/f-i70l-o©>";-;^i. 0 

[09] 01©l|g©*fC08©Ky7 h«S!SI4iffliitf 
[010] ASIU^BB. *^WiCj:S08©K'J7 

m^mmri,::^7' v ^Ryj#g*^-r 7 n - 5^ + - h t? 

[01 1 ] 01 0©A&CXBK:ffl^r«ffl3n-5fia 
[0 1 2 ] 0 1 ©cpCCiHjijn/c: I -o©;^ 'j 5. h^© 

[0 1 3 ] 0 1 ©*K:iBji*n/c 1 -3©f^#W&x V 
h*S©ll]teS^«:*JW?,jyc?rS3SrS 7 + - 
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1. Htle of Invention 

METHOD m APPARATUS FOR ANALYZING SPECTRAL DATA AND COMPUTER READABLE 

STORAGE MEDIUM FOR UTILIZATION TO ANALYZE SPECTRAL DATA 

2. Claims 

1. A »Gthod for analyzing spoctral data in a spectrometric 
instriment including a dispersion element and a detector 
receptive of dispersed light from the element, the detector 
having a plurality of detecting subarcays, with each subarray 
being at a different position on the detector; the method 
comprising steps of; 

acquiring first spectral data for a drift standard for selected 
subarray positions at a first time; 

comparing the first spectral data to a preassigned zero position 
for each selected subarray to obtain first offset data; 

acquiring second spectral data for a drift standard for selected 
subarray positions at a second time; 

comparing the second spectral data to the zero position for each 
selected subarray to obtain second offset data; and 

utilizing the first offset data and the second offset data to 
obtain a spectral shift for any subarray position at any selected 
time relative to the first time. 

2. The method of clain X wherein the step of utilizing comprises 
using the first offset data to obtain a first offset function 
defining an offset for any subarray position, using the second 
offset data to obtain a second offset function defining an offset 
for any subarray position, and utilizing the difference between 
the first offset function and the second offset function to 
obtain the spectral shift. 



(29) 



^mW- 1 0- 1 42054 



3. The nethod of claim a wherein the drift standard for the 
second epectral data is the drift standard Cor the first spectral 
data, and the selected subarray positions for the second spectral 
data are the selected subarray positions for the first spectral 
data. 

4. The method of claim 1 wherein the selected time is between 
the first tine and the second tine, and the spectral shift is 
obtained by interpolation between the first offset function and 
the second offset function. 

5. The method of claim i wherein the detector has a detector 
surface, and the detecting subarrays constitute a small portion 
of the detector surface. 

6. The method of claim 5 wherein the solected subarrays are 
substantially fewer in number than the plurality of detecting 
subarrays . 

7. The method of claim 1 further comprising acquiring test 
spectral data for a test sample at the selected time, and using 
the spectral shift to standardize the test spectral data to 
hypothetical instrument conditions. 

8. The method of claim 1 further comprising the steps of 
obtaining a base matrix model formed of base spectral data for at 
least one selected analyte, acquiring test .spectral data for a 
test sample at the selected time, using the spectral shift to 
shift the base spectral data so as to effect a shifted matrix 
model associated with instrument conditions at the selected time, 
and applying the shifted matrix model to the test spectral data 
so as to yield a parameter representing concentration of the 
selected analyte. 
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9. The method of claim 1 wherein, to obtain a base matrix model 
formed of base spectral data for at least one selected analyte, 
the matrix mod^l being for use with the test spectral data so as 
to yield a parameter representing concentration of the selected 
analyte, the method further comprises identifying the first time 
as an initial time/ the second time as a subsequent time, and the 
selected time as an interim time, obtaining preliminary spectral 
data for each selected analyte at the interim time, effecting the 
steps of acquiring, comparing and utilizing to obtain the 
spectral shift as an interim spectral shift for the interim time, 
and applying the interim spectral shift to the preliminary 
spectral data to effect the base spectral data for the base 
matrix model associated with instrument conditions at the initial 
time. 

10. The method of claim 9 wherein the test spectral data are 
acquired at predetermined increments in associated subarrays, the 
preliminary spectral data are obtained for a multiplicity of sub- 
increments smaller than the predetermined increments, and the 
shifted matrix model is applied with model spectral data culled 
for the predetermined increments from the preliminary spectral 
data. 

11. The method of claim 10 wherein the subarrays are each formed 
of a plurality of photosensitive pixels with a predetermined 
number of increments in each pixel, the increments being effected 
by slit-scanning with the instrument, 

12. The method of claim 11 wherein the multiplicity of sub- 
increments are effected by the slit-scanning, the predetermined 
number is generally non-integral, and the method further 
comprises steps of determining an axis shifting of spectral 
positions attributed to the predetermined number being non- 
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integral, and using the axis shifting to correct for the 
predeterriined number being non-integral. 

13, The method of claim 14 wherein the step of determining the 
axis shifting comprises further steps of: 

selecting an integral nominal number approximating the 
predetermined number, and a series of multiplicities of sub- 
increments for the nominal number and auxiliary numbers of sub- 
increments, the auxiliary numbers being smaller and larger than 
in the nominal total; 

acquiring spectral data for a selected spectral feature across a 
preselected pixel for each of the selected multiplicities to 
effect an associated original data series; 

shifting each data series in spectral position by one pixel to 
effect corresponding shifted data series; 

subtracting each original shifted data series from its 
corresponding shifted data series to effect a first set of 
differences for the smaller auxiliary numbers and a second set of 
differences for the larger auxiliary numbers; 

fitting the first set of differences to a first straight line and 
the second set of differences to a second straight line; and 

ascertaining an intersection point for the first straight line 
and the second straight line, whereby the intersection point has 
a displacement determinative of the axis shifting. 

14. The method of claim 12 wherein the step of using comprises 
determining a correction factor from the axis shifting, and 
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applying the correction factor to the base spectral data to 
effect corrected base spectral data for the base matrix model. 

15, The method of claim 13 wherein spectral data are related by 
a function 

= aa'Sfo + b(dYo/dXo) + e 
where Yq is derivative spectral data for sub-increnents having an 
increment size corresponding to the tentative total, Yj is 
increment spectral data for the sub-increments, is spectral 
position in the sub-increments , dYo/dX, is a derivative, a and b 
are parameters such that a parameter ratio b/a represents a 
tentative axis shift, and e is a fitting residual; the method 
comprising steps of: 

(a) estimating a tentative sub-increment sise corresponding to an 
estimated tentative total for each of the selected subarrays; 

(b) acquiring derivative spectral data, and further acquiring 
increment spectral data using the tentative increment size, for 
each of the selected subarrays; 

(c) fitting the derivative spectral data and the increment 
spectral data to the function to compute the parameter ratio for 
each of the selected subarrays; 

(d) fitting the increment size and the parameter ratio to a curve 
to ascertain deviation of the curve from a straight line; 

(e) using the deviation to estimate a corrected magnification 
corresponding to axis shifting; and 

if) repeating steps (b) through (e) until any deviation from a 
straight line in step <d) is less than a preselected limit, 
thereby effecting the integral total- 
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16. The method of claim 1 wherein the drift standard is a 
standard sample containing at least one analyte to effect a 
spectral peak in each of the selected subarray positions, 

17. The method of claim 1 wherein the drift standard is an 
optical element receptive of a light source to effect a series of 
regular secondary spectral peaks related to spectral positions 
within each of the selected subarray positions. 

18. The method of claim 17 wherein the instrument has a nominal 
calibration for spectral position versus spectral positions in 
the eubarrays, the instrument further includes a calibration 
source of a primary spectral peak having an identified spectral 
position, each secondary peak has an integer order number 
identified by correlation function to a peak spectral position in 
accordance with a correlation constant and a predetermined index 
of refraction of the interference element; and, to relate the 
secondary spectral peaks to spectral positions, the method 
further comprises: 

(a) acquiring primary spectral data for the primary peak in a 
first subarray position, and secondary spectral data for 
secondary peaks in the first subarray position and in other 
selected subarray positions; 

(b) estimating an initial correlation constant and a specified 
range thereof; 

(c) with the correlation function, the identified spectral 
position, and the initial correlation constant, calculating a 
tentative order number, and selecting a nearest integer order 
number to the tentative order number; 
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(d> with the correlation function, the idontificd spectral 
position and the nearest integer order number, calculating a 
corresponding first correlation constant to thereby effect a 
nuanber-constant pair consisting of the nearest integer order 
number and the first correlation constant; 

(e) recalculating a corresponding correlation constant vith a new 
order number constituting said nearest integer order nunber 
shifted by one to thereby effect a further number-constant pair 
consisting of the shifted order number and the corresponding 
correlation constant; 

(f) repeating step (e) with further integer order numbers shifted 
by additional ones until a full set of number-constant pairs is 
effected for the specified range of correlation constant; 

(g) in another selected subarray, identifying a secondary 
spectral peak of the secondary spectral data to a nominal peak 
spectral position determined by the nominal calibration; 

(h) with the correlation function, the nominal peak spectral 
position and each correlation constant of the set, computing 
further order numbers to effect additional number -const ant pairs; 

(i) designating all correlation constants in the additional pairs 
associated with order numbers that are non-integers, and deleting 
all number-constant pairs from the full set having the designated 
correlation constants^ thereby reducing the set of number- 
constant pairs; 

(j) repeating steps (g) , (h) and (i) until a single correlation 
constant in the pairs of the set remains to establish an 
effective correlation constant and corresponding remaining 
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integer order numbers for the calibration subarray and each 
selected subarray; and 

()c) with the correlation runction, the effective correlation 
constant and the remaining integer order numbers in the set, 
computing the spectral position for each selected secondary peak 
in each selected subarray. 

19. The method of claim 18 wherein the index of refraction is 
dependent on temperature and subarray position, the primary peak 
has a spectral position representative of temperature, and the 
method further comprises: 

while maintaining the interference element at a first 
temperature, effecting steps (a) through (k) with an estimated 
nominal index of refraction to ascertain a first primary peak 
position, a first effective correlation constant and associated 
order numbers; 

computing first values of the index of refraction with the 
correlation function using the first primary peak position, the 
first effective correlation constant and associated order 
numbers ; 

while maintaining the interference element at a second 
temperature, effecting steps (a) through (k) with the nominal 
index to ascertain a second primary peak position, a second 
effective correlation constant and associated order numbers; 

computing second values of the index of refraction with the 
correlation function using the second primary peak position, the 
second effective correlation constant and associated order 
numbers, the second values being dependent on the subarray 
positions; 
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computirg a value difference between each first value of the 
index and its correcponding second value for each corresponding 
subarray position, a position difference between the first 
primary pealc position and the second primary peaX position, and a 
difference ratio of each value difference to the position 
difference, and? 

storing the resulting difference ratios as a function of subarray 
position for subsequent use in computing the spectral position 
for each selected secondary peak in each selected subarray. 

20. The method of claim 19 further comprising: 

acquiring subsequent prinary spectral data defining an associated 
subsequent primary peak position at any selected time associated 
with a subsequent temperature of the interference element; 

computing a subsequent difference between the subsequent primary 
peak position and the first primary peak position, a 
multiplication product of the subsequent difference and each 
difference ratio to effect corrections in the index of 
refraction, and totals of the changes and the nominal index of 
refraction to effect a temperature corrected index of refraction 
for the selected subarray positions; and 

while maintaining the interference element at the subsequent 
temperature, effecting step (K) with the corrected index of 
refraction to compute the spectral position for each selected 
secondary peak in each selected subarray. 

21. A method for determining spectral position for a selected 
secondary peak for an optical interference element in a 
spectrometric instrument, the instrument including a dispersion 
element and a detector receptive of dispersed light from the 
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element, the detector having a plurality of detecting subarrays, 
each subarray being at a different position on the detector, the 
interference element being receptive of a light source to effect 
through the dispersion element and the detector a series of 
regular secondary spectral peaks related to spectral positions in 
the subarrays, the instrument having a nominal calibration for 
spectral position versus spectral positions in the subarrays, the 
instrument further including a calibration source of a primary 
spectral peak having an identified spectral position, each 
secondary peak having an integer order number identified by a 
correlation function to a peak spectral position in accordance 
with a correlation constant and a predetermined index of 
refraction of the Interference element; wherein, to relate the 
secondary spectral peaks to spectral positions, the method 
comprises: 

(a) acquiring primary spectral data for the primary peak in a 
first subarray position, and secondary spectral data for 
secondary peaks in the first subarray position and in other 
selected subarray positions; 

(b) estimating an initial correlation constant and a specified 
range thereof; 

(c) with the correlation function, the identified spectral 
position, and the initial correlation constant, calculating a 
tentative order number, and selecting a nearest integer order 
number to the tentative order number; 

(d) with the correlation function i the identified spectral 
position and the nearest integer order number, calculating a 
corresponding first correlation constant to thereby effect a 
number- const ant pair consisting of the nearest integer order 
number and the first correlation constant; 
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(e) recalculating a corresponding correlation constant with a new 
order number constituting said nearest integer order number 
shifted by one .to thereby effect a further number-constant pair 
consisting of the shifted order number and the corresponding 
correlation constant; 

(f) repeating step (e) with further integer order numbers shifted 
by additional ones until a set of number-constant pairs is 
effected for the specified range of correlation constant; 

(g) in another selected subarray, identifying a secondary 
spectral peak of the secondary spectral data to a nominal peak 
spectral position determined by the nominal calibration,- 

(h) with the correlation function, the nominal peak spectral 
position and each correlation constant of the set, computing 
further order numbers to effect additional number-constant pairs; 

(i) designating all correlation constants in the additional pairs 
associated with order numbers that are non- integers, and deleting 
all number-constant pairs from the full set having the designated 
correlation constants, thereby reducing the set of number- 
constant pairs; 

(j) repeating steps (g) , (h) and (i) until a single correlation 
constant in the pairs of the set remains to establish an 
effective correlation constant and corresponding remaining 
integer order numbers for the calibration subarray and each 
selected subarray; and 

(k) with the correlation function, the effective correlation 
constant and the remaining integer order numbers in the set, 
computing the spectral position for each selected secondary peak 
in each selected subarray. 
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22, The method of claim 21 wherein the index of refraction is 
dependent on temperature and subarray position, the primary peak 
has a spectral position representative of temperature, and the 
method further comprises: 

while maintaining the interference element at a first 
temperature, effecting steps (a) through (k) with an estimated 
nominal index of refraction to ascertain a first primary peak 
position, a first effective correlation constant and associated 
order numbers ; 

computing first values of the index of refraction with the 
correlation function using the first primary peak position, the 
first effective correlation constant and associated order 
numbers ; 

while maintaining the interference element at a second 
temperature, effecting steps (a) through (k) with the nominal 
index to ascertain a second primary peak posit ion« a second 
effective correlation constant and associated order numbers; 

computing second values of the index of refraction with the 
correlation function using the second primary peak position, the 
second effective correlation constant and associated order 
numbers, the second values being dependent on the subarray 
positions; 

computing a value difference between each first value of the 
index and its corresponding second value for each corresponding 
subarray position, a position difference between the first 
primary peak position and the second primary peak position, and a 
difference ratio of each value difference to the position 
difference; and 
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storing the resulting differencB ratios as a function of subarray 
position for subsequent use in computing the spectral position 
for each selected secondary peak in each selected subarray. 

23. The method of claim 22 further comprising: 

acquiring subsequent primary spectral data defining an associated 
subseq[uent primary peak position at any selected time associated 
with a subsequent temperature of the interference element; 

computing a subsequent difference between the subsequent primary 
peak position and the first primary peak position, a 
multiplication product of the subsequent difference and each 
difference ratio to effect corrections in the index of 
refraction, and totals of the changes and the nominal index of 
refraction to effect a temperature corrected index of refraction 
for the selected subarray positions; and 

while maintaining the Interference element at the subsequent 
temperature, effecting steps (a) through (k) with the corrected 
index of refraction to compute the spectral position for each 
selected secondary peak in each selected subarray. 

24. An apparatus for analyzing spectral data, including a 
spectrometric instrument having a dispersion element and a 
detector receptive of dispersed light fron the element, the 
detector having a plurality of detecting subarrays, with each 
subarray being at a different position on the detector; the 
apparatus comprising: 

means for acquiring first spectral data for a drift standard for 
selected subarray positions at a first time; 



< 
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means for comparing the first spectral data to a preassigned zero 
position for each selected subarray to obtain first offset data; 

neans for acquiring second spectral data for a drift standard for 
selected subarray positions at a second time; 

means for comparing the second spectral data to the zero position 
for each selected subarray to obtain second offset data; and 

means for utilizing the first offset data and the second offset 
data to obtain a spectral shift for any subarray position at any 
selected time relative to the first timfi- 

25. The apparatus of claim 24 wherein the means for utilizing 
comprises means for using the first offset data to obtain a first 
offset function defining an offset for any subarray position, 
means for using the second offset data to obtain a second offset 
function defining an offset for any subarray position, and means 
for utilizing the difference between the first offset function 
and the second offset function to obtain the spectral shift. 

26. The apparatus of claim 24 wherein the drift standard for the 
second spectral data are the drift standard for the first 
spectral data, and the selected subarray positions for the second 
spectral data are the selected subarray positions for the first 
spectral data. 

27. The apparatus of claim 24 wherein the selected time is 
between the first time and the second time, and the spectral 
shift is obtained by interpolation between the first offset 
function and the second offset function. 
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28. The apparatus of claim 24 wherein the detector has a 
detector surface, and the detecting subarrays constitute a small 
portion of the- detector surface. 

29. The apparatus of claim 28 wherein the selected subarrays are 
substantially fewer in number than the plurality of detecting 
subarrays. 

30. The apparatus of claia 24 further comprising means for 
acquiring test spectral data for a test sample at the selected 
time, and means for using the spectral shift to standardize the 
test spectral data to hypothetical instrument conditions. 

31. The apparatus of claim 24 further comprising a base matrix 
model formed of base spectral data for at least one selected 
analyte, means for accjuiring test spectral data for a tost sample 
at the selected time, means for using the spectral shift to shift 
the base spectral data so as to effect a shifted matrix model 
associated with instrument conditions at the selected time, and 
means for applying the shifted matrix model to the test spectral 
data so as to yield a parameter representing concentration of the 
selected analyte. 

32. The apparatus of claim 24 wherein, to obtain a base matrix 
model formed of base spectral data for at least one selected 
analyte, the matrix model being for use with the test spectral 
data so as to yield a parameter representing concentration of the 
selected analyte, the first time is an initial time, the second . 
tine is a subsequent time, and the selected time is an interim 
time, the apparatus further comprises means for obtaining 
preliminary spectral data for each selected analyte at the 
interim time, means for effecting the steps of acquiring, 
comparing and utilizing to obtain the spectral shift as an 
interim spectral shift for the interim time, and means for 
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applying the interin spectral shift to the preliminary spectral 
data to effect the base spectral data for the base matrix model 
associated with instrument conditions at the initial tine. 

33. The apparatus of claim 32 vherein the test spectral data are 
acquired at predetemined increments in associated subarrays, the 
preliminary spectral data are obtained for a multiplicity of sub- 
increments smaller than the predetermined increnents, and the 
shifted matrix model is applied with model spectral data culled 
for the predetermined increments from the preliminary spectral 
data. 

34. The apparatus of claim 33 wherein the subarrays are each 
formed of a plurality of photosensitive pixels with a 
predetermined number of increments in each pixel, the increments 
being effected by slit-scanning with the instrument. 

35. The apparatus of claim 34 wherein the multiplicity of sub- 
increments are effected by the slit- scanning, the predetemined 
number is generally non-integral, and the apparatus further 
comprises means for determining an axis shifting of spectral 
positions attributed to the predetermined number being non- 
integral, and means for using the axis shifting to correct for 
the predetermined number being non-integral. 

36. The apparatus of claim 33 wherein the means for determining 
the axis shifting comprises: 

means for acquiring spectral data for a selected spectral feature 
across a preselected pixel for each of a series of multiplicities 
of sub-increments, one xiultiplicity having a preselected integral 
nominal number of sub-increments approximating the predetermined 
number, and other multiplicities having auxiliary numbers of sub- 
increments, the auxiliary numbers being smaller and larger than 
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the nominal total, such spectral data affecting an associated 
original data series; 

neane for shifting «ach data series in spectral position by one 
pixel to effect corresponding shifted data series; 

means for subtracting each original shifted data series from its 
corresponding shifted data series to effect a first set of 
differences for the smaller auxiliary numbers and a second sot of 
differences for the larger auxiliary numbers; 

means for fitting the first set of differences to a first 
straight line and the second set of differences to a second 
straight line; and 

means for ascertaining an intersection point for the first 
straight line and the second straight line, whereby the 
intersection point has a displacement determinative of the axis 
shifting. 

37. The apparatus of claim 35 wherein the means for using 
comprises means for determining a correction factor from the axis 
shifting, and means for applying the correction factor to the 
base spectral data to effect corrected base spectral data for the 
base matrix model. 

38. The apparatus of claim 35 wherein spectral data are related 
by a function 

Y| - ajYo + b<dV./dXo) + e 
where Y, is derivative spectral data for sub- increments having an 
increment sise corresponding to the tentative total, Y, is 
increment spectral data for the cub-incremants , Ij is spectral 
position in the sub-increments, dY,/dXa is a derivative, a and b 
are parameters such that a parameter ratio b/a represents a 
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tentative axis shift, and o is a fitting residual; and the 
apparatus further comprises the stored function, and further 
comprises: 

(a) means for estimating a tentative sxib- increment size 
corresponding to an estimated tentative total for each of the 
selected subarrays; 

(b) means for acquiring derivative spectral data, and further 
acquiring increment spectral data using the tentative increment 
size, for each of the selected subarrays; 

(c) means for fitting the derivative spectral data and the 
increnent spectral data to the function to compute the parameter 
ratio for each of the selected subarrays; 

(d) means for fitting the increment size and the parameter ratio 
to a curve to ascertain deviation of the curve from a straight 
line; 

(e) means for using the deviation to estimate a corrected 
magnification corresponding to axis shifting; and 

(f) means for repetitively applying the means (b) through (e) 
until any deviation from a straiglit line of means (d) is less 
than a preselected limit, thereby effecting the integral total, 

39. The apparatus of claim 24 wherein the drift standard is a 
standard sample containing at least one analyte to effect a 
spectral peak in each of the selectad subarray positions, 

40. The apparatus of claim 24 wherein the drift standard is an 
optical element receptive of a light source to effect a series of 
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regular secondary spectral peaks related to spectral positions 
within each of the selected subarray positions. 

41, The apparatus of claim 4 0 wherein the instrument has a 
nominal calibration for spectral position versus spectral 
positions in the subarrays, the instrunent further includes a 
calibration source of a primary spectral peak having an 
identified spectral position, each secondary peak has an integer 
order number identified by correlation function to a peak 
spectral position in accordance with a correlation constant and a 
predetermined index of refraction of the interference element 
and, to relate the secondary spectral peaks to spectral 
positions, the apparatus further comprises: 

(a) means for acquiring primary spectral data for the primary 
peak in a first subarray position, and secondary spectral data 
for secondary peaks in the first subarray position and in other 
selected subarray positions; 

(b) means for calculating a tentative order number with the 
correlation function, the identified spectral position, a pre- 
estinated initial correlation constant and a specified range 
thereof ; 

(c) means for selecting a nearest integer order number to the 
tentative order number; 

(d) means for calculating a corresponding first correlation 
constant with the correlation function, the identified spectral 
position and the nearest integer order number, so as to thereby 
effect a number-constant pair consisting of the nearest integer 
order number and the first correlation constant; 
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<e) means for recalculating a corresponding correlation constant 
with. a new order number constituting said nearest integer order 
number shifted. by one to thereby effect a further nuinber-constant 
pair consisting of the shifted order number and the corresponding 
correlation constant; 

(f ) means for repetitively applying the means (e) with further 
integer order numbers shifted by additional ones until a set of 
number-constant pairs is effected for the specified range of 
correlation constant; 

(g) means for identifying, in another selected subarray, a 
secondary spectral peak of the secondary spectral data to a 
nominal peak spectral position determined by the nominal 
calibration; 

(h) means for computing further order numbers to effect 
additional number-constant pairs, said means utilizing the 
correlation function, the nominal peak spectral position and each 
correlation constant of the set; 

(i) means for designating all correlation constants in the 
additional pairs associated with order numbers that are non- 
integers, and for deleting all number-constant pairs from the 
full set having the designated correlation constants, thereby 
reducing the set of number-constant pairs; 

(j) means for repetitively applying the means (g) , (h) and (i) 
until a single correlation constant in the pairs of the set 
remains to establish an effective correlation constant and 
corresponding remaining integer order numbers for the calibration 
subarray and each selected subarray; and 
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(X) means for computing the Bpactral position for each selected 
secondary peak in each selected subarray, said means utilizing 
the correlation, function, the effective correlation constant and 
the remaining integer order numtoers in the set. 

42. The apparatus of claim 41 wherein the index is dependent on 
temperature and subarray position, the primary peak has a 
spectral position representative of temperature, and the 
apparatus further comprises: 

means for repetitively applying the iieans (a) through (K) with a 
predetermined nominal index of refraction, while maintaining the 
interference element at a first temperature, so as to ascertain 
first primary spectral data defining an associated first primary 
peak position, a first effective correlation constant and 
associated order numbers ; 

means for computing first values of the index of refraction with 
the correlation function using the first primary peak position, 
the first effective correlation constant and associated order 
numbers ; 

weans for repetitively applying the means (a) through (k) with 
the nominal index, while maintaining the interference element at 
a second temperature, so as to ascertain second primary spectral 
data defining an associated second primary peak position, a 
second effective correlation constant and associated order 
numbers ; 

means for computing second values of the index of refraction with 
the correlation function using the second primary peak position, 
the second effective correlation constant and associated order 
numbers, the second values being dependent on the subarray 
positions; 
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means for computing a value difference between each first value 
of the index and its corresponding second value for. each 
corresponding subarray position, a position difference l>etween 
the first primary pea)c spectral position and the second primary 
peak spectral position, and a difference ratio of each value 
difference to the position difference; and 

means for storing the resulting difference ratios as a function 
of subarray position for subsequent use in computing the spectral 
position for each selected secondary peak in each selected 
subarray. 

43. The apparatus of claiia 42 further comprising: 

means for acquiring subsequent primary spectral data defining an 
associated subsequent primary peak position at any selected time; 

means for computing a subsequent difference between the 
subsequent primary peak position and the first primary peak 
position, a multiplication product of the sxibsequent difference 
and each difference ratio to effect corrections in the index of 
refraction, and totals of the changes and the nominal index of 
refraction, so as to effect a temperature corrected index of 
refraction for the selected subarray positions; and 

means for repetitively applying the means (a) through <k) with 
the nominal index, while maintaining the interference element at 
the subsequent temperature, so as to compute the spectral 
position for each selected secondary peak in each selected 
subarray. 

44, An apparatus for determining spectral position for a 
selected secondary peak for an optical interference element in a 
spectroaetric instrument, the instrument including a dispersion 
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element and a detector receptive of dispersad light from the 
element, the detector having a plurality of detecting subarrays, 
each subarray being at a different position on the detector, the 
interference element being receptive of a light source to effect 
through the dispersion element and the detector a series of 
regular secondary spectral peaks related to spectral positions in 
the subarrays, the instruaent having a noainal calibration for 
spectral position versus spectral positions in the subarrays, the 
instrument further including a calibration source of a prltaary 
spectral peaX having an identified spectral position, each 
secondary peak having an integer order number identified by a 
correlation function to a peak spectral position in accordance 
with a correlation constant and a predetermined index of 
refraction of the interference eleioent; wherein, to relate the 
secondary spectral peaks to spectral positions, the apparatus 
comprises: 

(a) means for acquiring primary spectral data for the primary 
peak in a first subarray position, and secondary spectral data 
for secondary peaks in the first subarray position and in other 
selected subarray positions; 

(b) means for calculating a tentative order number with the 
correlation function, the identified spectral position, a pre- 
estimated initial correlation constant and a specified range 
thereof; 

(c) means for selecting a nearest integer order number to the 
tentative order number; 

(d) means for calculating a corresponding first correlation 
constant with the correlation function, a predetermined index of 
refraction, the identified spectral position and the nearest 
integer order number, so as to thereby effect a number-constant 
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pair consisting of the nearest integar order nmaber and the first 
correlation constant; 

(e) means for recalculating a corresponding correlation constant 
with a new order nuinber constituting said nearest integer order 
number shifted by one to thereby effect a further number- constant 
pair consisting of the shifted order number and the corresponding 
correlation constant; 

(f ) means for repetitively applying the neanc (a) with further 
integer order numbers shifted by additional ones until a set of 
number- const ant pairs is effected for the specified range of 
correlation constant; 

(g) means for identifying, in another selected subarray, a 
secondary spectral peak of the secondary spectral data to a 
nominal peak spectral position determined by the nominal 
calibration; 

(h) means for computing further order numbers to effect 
additional number-constant pairs, said means utilising the 
correlation function, the nominal peak spectral position and each 
correlation constant of the set; 

(i) means for designating all correlation constants in the 
additional pairs associated with order numbers that are non- 
integers, and for deleting all number -constant pairs from the 
full set having the designated correlation constants, thereby 
reducing the set of number-constant pairs; 

(j) means for repetitively applying the means (g) , (h) and (i) 
until a single correlation constant in the pairs of the set 
remains to establish an effective correlation constant and 
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corresponding remaining integer order numbers for the calibration 
subarray and each ealected subarray; and 

(k) raeans for computing the spectral position for each selected 
secondary peaX in each selected subarray, said Jueans utilizing 
the correlation function, the effective correlation constant and 
the remaining integer order numbers in the set. 

45. a:he apparatus of claim 44 wherein the index is dependent on 
temperature and subarray position, the primary peak has a 
spectral position representative of temperature, and the 
apparatus further comprises: 

means for repetitively applying the means (a) through (k) with a 
predetermined noninal index of refraction, while maintaining the 
interference element at a first temperature, so as to ascertain 
first primary spectral data defining an associated first primary 
peak position, a first effective correlation constant and 
associated order numbers ; 

means for computing first values of the index of refraction with 
the correlation function using the first primary peak position, 
the first effective correlation constant and associated order 
numbers ; 

means for repetitively applying the means (a) through (k) with 
the nominal index, while maintaining the interference element at 
a second temperature, so as to ascertain second primary spectral 
data defining an associated second primary peak position, a 
second effective correlation constant and associated order 
numbers; 

means for computing second values of the index of refraction with 
the correlation function using the second primary peak position, 
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the second effective correlation constant and associated order 
numbers, the second values being- dependent on the subarray 
positions ; 

means for computing a value difference between each first value 
of the index and its corresponding second value for each 
corresponding subarray position, a position difference between 
the first primary peak spectral position and the second primary 
peak spectral position, and a difference ratio of each value 
difference to the position difference; and 

means for storing the resulting difference ratios as a function 
of subarray position for subsequent use in computing the spectral 
position for each selected secondary peak in each selected 
subarray. 

46. The apparatus of claim 4 5 further comprising: 

means for acquiring subsequent primary spectral data defining an 
associated subsequent primary peak position at any selected time; 

means for computing a subsequent difference between the 
subsequent primary peak position and the first primary peak 
position, a multiplication product of the subsequent difference 
and each difference ratio to effect corrections in the index of 
refraction, and totals of the changes and the nominal index of 
refraction, so as to effect a temperature corrected index of 
refraction for the selected subarray positions; and 

means for repetitively applying the means (a) through (k) with 
the nominal index, while maintaining the interference element at 
the subsequent temperature, so as to compute the spectral 
position for each selected secondary peak in each selected 
subarray. 
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47- A computer readable storage medium for utilization to 
analyze spectral data for a sanplo in a spectroBfttric instrunent 
that includes a dispersion elesient and a detsctor receptiva of 
dispersed light from tlx* element, the detector having a plurality 
of detecting cubarrays with each subarray being at a different 
position on thft detector, the InstruBient further including means 
for acquiring first spectral data for a drift standard for 
selected subarray positions at a first time, means for acquiring 
second spectral data for a drift standard for selected subarray 
positions at a second time, and computing means receptive of the 
spectral data for computing corresponding spectral information 
representative of the sample, the storage medium having data code 
and program code embodied therein so as to be readable by the 
computing means, wherein the data code comprises a preassigned 
zero. position for each selected subarray, and the program code 
comprises means for comparing the first spectral data to the 
preassigned aero position for each selected subarray to obtain 
first offset data, means for comparing the second spectral data 
to the zexQ position for each selected subarray to obtain second 
offset data, and means for utilizing the first offset data and 
the second offset data to obtain a spectral shift for any 
subarray position at any selected time relative to the first 
time. 

48. The storage medium of claim 47 vherein the means for 
utilizing comprises means for using the first offset data to 
obtain a first offset function defining an offset for any 
subarray position, means for using the second offset data to 
obtain a second offset function defining an offset for any 
subarray position, and means for utilizing the difference batueen 
the first offset function and the second offset function to 
obtain the spectral shift. 
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49. The storage nediua of claim 47 wherein th» fialacted time is 
between tlxe first tljae and tlx© caoond time, and the spectral 
shift is obtained by interpolation between the first offset 
function and the second offset function. 

50. The storage mediun of claim 47 wherein the instrument 
further includes means for acquiring test spectral data for a 
test sariple at the selected time, and the program code further 
comprises means for using the spectral shift to standardize the- 
test spectral data to hypothetical instrument conditions. 

51. The storage medium of claim 47 wherein the instrument 
further includes means for acquiring test spectral data for a 
test sample at. the selected time, the data code further comprises 
a base matrix model formed of base spectral data for at least one 
selected* analyte, and the program code further comprises means 
for using the spectral shift to shift the base spectral data so 
as to effect a shifted matrix model associated with instrument 
conditions at the selected time, and means for applying the 
shifted matrix model to the test spectral data so as to yield a 
parameter representing concentration of the selected analyte. 

52- The storage medium of claim 47 wherein the instrument has a 
nominal calibration for spectral position versus spectral 
positions in the subarrays, the instrument further includes a 
calibration source of a primary spactral peak having an 
identified spectral position, the drift standard is an optical 
element receptive of a light source to effect a series of regular 
secondary spectral peaJcs related to spectral positions within 
each of the selected subarray positions, each secondary peak has 
an integer order number identified by correlation function to a 
peaX spectral position in accordance with a correlation constant 
and a predetairmined index of refraction of the interference 
element, and the Instrument further includes: 
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(a) moans for acijuiring primary spectral data tor the primary 
peak in a first subatxay position, and secondary spectral data 
tor secondary peaks in the first subarray position and In other 
selected subarray positions; and 

the progxBS codo further comprises : 

(b) means for calculating a tentative order number with the 
correlation function, the identified spectral position, a pre- • 
estiiaatea initial correlation constant and a specified range 
thereof; 

(c) means for selecting a nearest integer order number to the 
tentative order nximber; 

(d) aeans for calculating a corresponding first correlation 
constant with the correlation function, the identified spectral 
position and the nearest integer order nunber, so as to thereby 
effect a nuBber-constant pair consisting of the nearest integer 
order number and the first correlation constant; 

(e) Means for recalculating a corresponding correaation constant 
with a new order nuaber constituting said nearest Integer order 
number shifted by one to thereby effect a further number-constant 
pair consisting of the shifted order number and the corresponding 
correlation constant; 

(f ) means for repetitively applying the means (e) with further 
integer order numbers shifted by additional ones until a set of 
number-constant pairs is effected for the specified range of 
correlation constant; 

(g) means for identifying, in another selected subarray, a 
secondary spectral peaJc of the secondary spectral data to a 
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noninal peaX spectral position determined by the nominal 
calibration; 

(h) means for computing further order niuabers to effect 
additional number-constant pairs, said neans utilizing the 
correlation function, the nominal peak spectral position and each 
correlation constant of the set; 

(i] maans for designating all corralation constants in the 
additional pairs associated with order numbers that are non- 
integers, and for deleting all number-constant pairs froa the 
full set having the designated correlation constants, thereby 
reducing the set of number-constant pairs; 

(j) means for repetitively applying the means (g) , W and (i) 
until a single correlation constant in the pairs of the set 
remains to establish an effective correlation constant and 
corresponding remaining integer order numbers for the calibration 
subarray and each selected subarray; and 

(3c) means for computing the spectral position, for each selected 
secondary peak in each selected subarray, said means utilizing 
the correlation function, the effective correlation constant and 
the remaining integer order numbers In the set- 

53. x computer readable storage medium for utilization to 
determine spectral position for a selected secondary p«a)c for an 
optical interference element in a spectrometric instrument that 
includes a dispersion element and a detector receptive of 
dispersed light from the element, the detector having a plurality 
of detecting subarrays, each subarray buing at a different 
position on the detector, the interference element being 
receptive of a light source to effect through the dispersion 
element and the detector a series of regular secondary spectral 
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peaJcs related to spectral positions in tha subarrays, the 
instrument further including a calibration sourca of a prinary 
spectral peak having an integer order number identified by a 
correlation function to a peak spectral position in accordance 
with a correlation constant and a predetennined index of 
refraction of tue interference element, and the instrument 
further including: 

(a) means for accjuiring primary spectral data for the primary 
peak in a first fiubarray position, and secondary spectral data 
for secondary peaks in the first subarray position and in other 
selected subarray positions; 

the storage nedium having data cod© and program code embodied 
therein so as to be readable by tha computing means, wherein the 
data code coiapricoc a nominal calibration of the instrument for 
spectral position versus spectral positions in the subarrays, tha 
correlation function, and the predetermined index of refraction; 
and the program code comprises; 

(b) means for calculating a tentative order nuniber with the 
correlation function, the identified spectral position, a pre- 
estimated initial correlation constant and a specified range 
thereof; 

(c) means for selecting a nearest integer order number to the 
tentative order number; 

(d) means for calculating a corresponding first correlation 
constant with the correlation function, a predetermined irtdex of 
refraction, the identified spectral position and the nearest 
integer order number, so as to thereby effect a number-constant 
pair consisting of the nearest integer order number and the first 
correlation constant; 
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(e) means for recalculatiJig a coxre spoil ding correlation constant 
with a nev order nuiaber constituting said nearest integer order 
nuabar shifted by one to thereby effect a further nuiaber-constant 
pair consisting of the shifted order number and the correcponding 
correlation constant; 

(f } means for repetitively applying the means (a) vith further 
integer order numbers shifted by additional ones until a set of 
number-constant pairs is effected for the specified range of 
correlation constant; 

(g) means for identifying, in another eeleeted eubarray, a 
secondary spectral pea}c of the secondary spectral data ^o a 
nominal peaH spectral position determined by the nominal 
calibration; 

(h) means for computing further order numbers to effect 
additional number-constant pairs, said means utilizing the 
correlation function, the nominal peak spectral position and each 
correlation constant of the set; 

(i) moans for designating all correlation constants in the 
additional pairs associated with order numbers that are non- 
integers, and for deleting all number- constant pairs from the 
full set having the designated correlation constants, thereby 
reducing the set of number- constant pairs; 

(j) means for repetitively applying the means (g) , (h) and (i) 
until a single correlation constant in the pairs of the set 
remains to establish an effective correlation constant and 
corresponding remaining integer order numbers for the calibration 
subarray and each selected subarray; and 



()c) means for computing the spectral position for each selected 
secondary peak in each selected subarray, said means utilizing 
the correlation function, the effective correlation constant and 
the remaining integer order numbers in the eet. 
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3. Detailed Explanation of the Invention 
This invention relates to spectroaetric instruments and particularly to 
method and apparatus for analyiing spectral data and conputer readable 
storage mediuni for utilization to analyze spectral data for a sample in 
such instruments. 

BACXGROUKD 

Spectrometric instruments include a disper&ion eleaant, such as a 
diffraction grating, and a detector eystem. Modem instruments 
include a computer that is receptive of spectral data from the 
detector to analyze and compare spectra. With improvements in 
optics, detectors and computerization, there has evolved an 
ability to perform very precise meastirenents. 

One type of instriuaent utilizes an inductively coupled plasma 
(ZCF) with sample injection to effect spectral lines of atomic 
elements. A spectrophotometer used in conjunction with ICP is a 
crossed grating type that produces a two dimensional array of 
spectral lines. A detector for the array has segmented subarrays 
of photosensitive pixels located strategically only at expected 
locations of the spectral lines. A solid state device such as a 
charge coupled device (CCD) with the photosensitive pixels in the 
subarrays is used. A crossed grating spectrophotometer 
incorporating such a detector is disclosed In U.S. patent Ko. 
4 , 820, 048 (Barnard) . 

With the evolving requirements for precision, variations among 
instruments, and drift in each instrument (e.g. from temperature 
and pressure variations) have become more of a problem. U.S. 
patent Mo. 5,303,165 (Ganz et al, "Ganz patent") discloses 
standardization of instruments by transforming spectral data with 
a transformation filter derived from a specified line profile 
common to the instruments. Spectral lines as displayed by a 
spectrometric instrument actually have a finite width and 
profile, and tihe standardization corrects for variations in 
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instrument profile. Such standardization is distinguished from 
calibration associated with quantification. Determination of 
conpositional quantities of a sample is carried out separate!/ or 
in conjunction with standardization. Wavelength calibration may 
be associated vith standardization, but actual wavelengths of 
measured spectra are not necessarily needed. The technique of 
the Ganz patent is quite precise and useful, particularly with 
continuous array detectors, but suffers from a requirement for 
substantial amounts of spectral data collection and associated 
lengthy computations. The Ganz patent also discloses tho use of 
a source of regular fringe peaks for wavelength calibration, but 
such use is not suited for a segmented subarray detector. 

The presence and quantity of components in a sample may be 
determined with computer computations by application of 
calibration models to spectral data, the models being derived 
from spectra of known quantities of sample analytes (individual 
atomic elements) . An archive of model data is stored in computer 
memory for application to sample data for essentially automatic 
determination of components and their quantities in the sample. 
One approach is disclosed in U.S. patent No. 5,308,982 (Ivaldi et 
al) which incorporates a derivative of sample spectral data into 
the matrix model to compensate for spectral drift. This is a 
standardization that requires spectral data to be acquired in 
relatively small spectral increments to achieve sufficient 
representation ot the derivative in the model, wavelength 
increments of spectral data ordinarily is limited by pixel size. 
Smaller increments are achieved by slit scanning in which the 
inlet slit to the spectrometer is imaged on a pixel. Varying the 
lateral position of the slit in small steps effectively moves a 
spectrum across the pixels to obtain spectral data in smaller 
increments. Although utilized for collecting archive data, it is 
preferable that slit scanning be avoided to speed up ordinary 
data acquisition. 
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SUMMARY 

Objects of the invention are to provide a novel method and a 
novel means tor analyzing spectral data in a spec trome trie 
instrument having a plurality of detecting subarrays. other 
objects are to provide such a wethod and a means for correcting 
for instrument variations including drift. Other objects are to 
provide such a method and a means for adjusting model calibration 
data to instrument conditions currently with sample data 
acquisition. Further objects are to provide a novel method and a 
novftl means for constructing a calibration model for a 
spectrometrlc instrument. Further objects are to provide a novel 
method and a novel means for adjusting slit scanning in a 
spectrometric instrument. Additional objects are to provide a 
novel method and a novel means for utilizing a source of regular 
fringe peaks for wavelength calibration with a segmented subarray 
photodetactor in a spectrometrlc instrument. 
Further objects are to provide a conputer readable storage .ediu. (or 
utilization to analyte spectral data (or a sample in a spectrometric 
instrument. 



The foregoing and other objects are achieved, at least in part by 
analyzing spectral data in a spectrometric instrunent that 
includes a dispersion element and a detector receptive of 
dispersed light from the element, the detector having a plurality 
of detecting subarrays, with each subarray being at a different 
position on the detector. Advantageously the detecting subarrays 
constitute a small portion of the detector surface. The selected 
subarrays may be substantially fewer in number than the plurality 
of detecting subarrays. 

First spectral data are acquired for a drift standard for 
selected subarray positions at a first time, and compared to a 
preasslgned zero position for each selected subarray to obtain 
first offset data. Second spectral data are acquired for the 
drift standard for the selected subarray positions at a second 
time, and compared to the rero position for each selected 
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subarray to obtain second offset data. The first offset data and 
the second offset data are utilized to obtain a spectraJ shift 
for any subarray position at any selected time relative to the 
first time by interpolation or extrapolation. Advantageously, 
the first offset data are used to obtain a first offset function 
defining an offset for any subarray position, the second offset 
data are used to obtain a second offset function defining an 
offset for any subarray position, and the difference between the 
first offset function and the second offset function is utilized 
to obtain the spectral shift. 

Although the shift may be applied to test data, advantageously 
the shift is applied to a matrix model used for converting the 
data to compositional information. For this, o base matrix model 
is formed of base spectral data for at least one selected 
analyte. Test spectral data are acquired for a test sample at 
the selected time, using the spectral shift to shift the base 
spectral data to the model spectral data for the selected time to 
effect a shifted matrix model associated with instrument 
conditions at the selected time. The shifted matrix model is 
applied to the test spectral data so as to yield a parameter 
representing concentration of the selected analyte. 

Archive data for the matrix model preferably is obtained in the 
foregoing manner to account for spectral shift during data 
acquisition. More preferably the data are acquired with slit 
scanning in the instrument to achieve sub- increments smaller than 
pixel size in the detector, utilizing a procedure to assure that 
there Is known spacing between scanning steps across one pixel. 

The drift standard may be a chemical standard of selected 
analytes, or an optical element such as an interference element, 
producing regular secondary peaks (fringes) related to spectral 
positions in each subarray. To utilize an interference element, 
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the instrument further includes a calibration source of a primary 
spectral peak having an identified spectral position. Each 
secondary peak has an integer order number identified by a 
correlation function to a peak spectral position in accordance 
with a single correlation constant (related to thickness of the 
element) . To relate the secondary peaks to spectral positions, 
primary spectral data are acquired for the primary peak in a 
first subarray position, and secondary spectral data are acquired 
for secondary peaks in the first subarray position and in other 
selected subarray positions. An initial correlation constant 
with a specified range thereof are estimated. A tentative order 
number is calculated with the correlation function, the 
identified spectral position, and the initial correlation 
constant; and a nearest integer order number to the tentative 
order number is selected. A corresponding first correlation 
constant is recalculated with the nearest integer order number, 
so as to thereby effect a first number-constant pair. Another 
correlation constant is calculated with a new order number 
constituting said nearest integer order number shifted by one to 
effect another number-constant pair. The calculations are 
repeated with further integer order numbers shifted by additional 
ones until a set of number-constant pairs is effected for the 
specified range of correlation constant, for the first subarray. 

In another selected subarray , a secondary peak is identified to a 
nominal peak position determined by the nominal calibration. 
With the correlation function, the nominal peak position and each 
correlation constant of the pair set, further order numbers are 
computed to effect additional number-constant pairs. All 
correlation constants in the additional pairs associated with 
order numbers that are non-integers are designated, and all such 
corresponding pairs are deleted from the full set. The foregoing 
three steps are repeated until a single correlation constant in 
the pairs of the set remains to establish an effective 
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correlation constant and corresponding remaining integer order 
nimbers. With the correlation function, the affectiv* 
correlation constant and the remaining integer order numbers in 
the set, the spectral position is coraputed for each selected 
secondary peak in each selected subarray. It temperature varies, 
a procedure should be used to adjust the Index of refraction of 
the interference element used in the correlation function. 

DETAILED OESCRIPTIOK 

The invention is utilized with a conventional or other desired 
spectronetric instrument 10 (JIG* l) for spectral analysis in the 
general range of infrared, visible and ultraviolet radiation, 
where spectral features are in narrow banda or spectral lines 
received on an array photodetector 12 from ICP light source lA 
via an entrance lens is, a slit i7, another focusing lens 19 and 
a dispersion element is such as a grating or prism. The 
instrument utilizes, for example, an Inductively coupled plasma 
18 with sample injection 2 0 for the light source to effect atomic 
or ionic emission lines, such as a Perkin-Elmer Optima 3000 
instrument. Alternatively, the Instrument may be adapted for 
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molecular spectroscopy such as the type described in the 
aforementioned U.S. patent Ho. 5,303,165 (Ganz) , but with 
substitution of. the present invention for the standardization 
described therein. 

The detector is (FIG. 2) in such an instrument typically is a 
self scanned photodiode array, a charge injection device or 
preferably a charge-coupled device (CCD) with photosensitive 
pixels 22 distributed over a two dinensional surface 24 as 
segmented subarrays 26 of pixels at selected positions of 
expected spectral lines as disclosed in the aforementioned U.S. 
patent No, 4,820,048 (Barnard) of the present assignee, 
incorporated in Its entirety herein by reference. The spectral 
lines in this case are effected by crossed dispersion elements 
such as a pair of gratings or a grating and a prism. In other 
applications the pixels in an array cover the entire detecting 
surface of a linear or two dimensional detector. 

Reference to "subarray" herein and in the claims means either 
selected groupings of pixels otherwise distributed over an entire 
linear or tvo-dinensional surface, or segmented subarrays of 
pixels at selected positions. The term "pixel position" refers 
conventionally to spectral position in the finite increments of 
an array photodetector . A segmented array detector may contain, 
for example, 7000 pixels in 24 5 subarrays with each subarray 
containing 10 to 40 pixels, generally about 16 pixels. Spectral 
data 21 (PIG. 3) for a spectral feature 2 3 obtained from the 
detector via a signal line 28 (pig. i) and stored in a memory 
section 30 of the computer 32 of the instrument is representative 
of photo energy received by each pixel 22 in a subarray 26 during 
a cycle period. An instrument often is calibrated for wavelength 
(or wavenumber) against pixel positions; however, as explained 
below, this is not necessary in some aspects of the present 
invention where the instrument provides compositional information 
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without need for wavelength data. The invention is applied in 
conjunction with the computer to compensate for instrument drift 
with time wherein a spectral feature shifts at least fractionally 
over the pixels (or in finite increiaents of spectral positions). 
(Unless otherwise stated, as used herein and in the claims the 
terms "wavelength" and "wavenumber" mean either actual wavelength 
or wavenumber or frequency, or spectral position (e.g. in pixel 
units) in a detector array or subarray.) 

The computer 32 generally is conventional such as a Digital 
Equipment Corporation model 5100 incorporated into the instrument 
by the manufacturer thereof. The computer should include a 
central processing unit (CPU) 32 with appropriate analog/digital 
converters (in and/or out as required) , memory sections 30 that 
may include a disX, a keyboard for operator input, and a monitor 
36 and/or a printer for display of the desired output such as 
sample composition. Programming is conventional such as with 
pc*+" generally incorporated into the computer by the 
manufacturer of the computer or instrument. Adaptations of the 
programming for the present invention will be readily recognized 
and achieved by those skilled in the art. 

A drift standard, in the present embodiment being a mixed-element 
standard sample of known, fixed elemental composition, is 
introduced into the instrument, e.g. injected 20 into the ICP 14. 
It is not necessary to know the quantitative amounts of the 
compositions of the elements in this standard, other than that 
they provide an adequate ratio of signal to noise. The atomic 
elements for the composition ore chosen for their distributed 
regions of spectral activity on the detector. An example is 
shown in Pis. 4, which shows a selected subset of the subarray 
positions 3 8 targeted by the proposed drift standard. Positions 
on the detector are defined by x and y coordinates of each 
selected subarray. A suitable standard may be formed, for 
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example, from 8 common elements yielding 17 spectral lines. The 
number of selected subarrays in the subset may be substantially 
fewer than the total number of subarrays (e.g. 17 in the subset 
compared with 245 total) . 

With reference to a flow chart (FIG. 5) , each subarray is 
assigned a reference point or zero position 48. The zero 
positions may be associated with hypothetical data which 
generally may be any arbitrarily selected spectral data for each 
pixel. Such data may be that acquired at the first time, or a 
theoretical or other good estimate of a design position and shape 
for a spectral feature at each array. Advantageously, however, 
the zero position simply is an arbitrary position of the subarray 
such as the center of the subarray. 

A first set of spectral data 40 (also yic, 6) is acquired 42 for 
the subset of selected subarray positions at a first time 43 for 
the drift standard 44. These data define pea}c locations 45 which 
are compared 46, generally by subtraction, with the zero 
positions 48 to obtain first offset data 50. These data 
generally will depend on position on the detector, and thereby 
are a function of the x-y coordinates 38 (PIG. 4) of the 
subarrays. The position coordinates nay be an average for the 
subarray or for each pixel, the former generally being 
sufficient. A mathematical offset function 52 ("first offset 
function") is determined (computed) 54 from the first offset data 
such that a "pixel offset" o (real line position relative to the 
zero position) is defined for all of the subarrays as a function 
of any x and y position on the detector. The function may be in 
the form of an equation, e.g. offset 0=ax+by+e where the 
parameters a, b and c define the function. The change in this 
mathematical function is used to predict the degree of spectral 
drift in pixel space (e.g. to as small as 0.01 pixel drift) as a 
function of any x and y position on the detector. Thus pixel 
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drift is defined as the change in th© pixel offset calculated as 
a function of x and y. Use of the offset equation simplifies 
standardization, process by enabling drift for all subarrays to be 
determined, although only a subset is measured. A higher order 
equation may be used, but this generally should not be necessary. 

To detect drift, a second set of spectral data 56 is .acquired 42 
for a drift standard 44 at a second, later time fiO. The second 
standard and selected array positions 38 preferably (but not 
necessarily) are the sane as the first. These data are compared 
4£ similarly with the zexo position 4B to obtain second offset 
data 62. The latter is used to obtain 54 a second offset 
function 64 defining an offset for any array position at the 
second time. The difference 66 between the first and second 
offset functions 52, 64 provides a difference function 69 
representing the spectral shift as a function of x-y position 38 
during the interval. By interpolation 70 or extrapolation a 
spectral shift 72 nay be obtained for any selected (third) time. 

Spectral data 76 for a test sample 78 (typically unknown) 
acquired at a selected time ("time 3") ao may be standardized 82 
to any base time such as the first time 43, by application 74 of 
the spectral shift 72. Tor better accuracy the selected time 
should be between the first and second tines for interpolation. 
This type of standardization is an option as shown by the broken 
lines in FIG. 5. The model (explained below, but unshifted) may 
be applied to the standardized data. 

The standardisation preferably is utilized in conjunction with 
application of a matrix model of calibration data applied to 
sample data. In such application, a conventional matrix model is 
obtained initially in a manner taught in the aforementioned U.S. 
patent No. 5,308,982 (Ivaldi et al, "Ivaldi patent"), 
incorporated herein by reference, and an article "Advantages of 
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Coupling Multivariate Data Reduction Techniques with Inductively 
Coupled Plasma Optical Emission Spectra", by J,C. Ivaldi and T,W. 
Barnard, Spectrochinica Acta, 48B, 1265 (1993). 

The Hiatrix is formed of archival spectral data representing 
accurate, base data for Xnown concentrations of one or a 
plurality of analytes expected to be found in test saaples. The 
matrix is formed of columns of spectral data, each column being 
for one analyte. The matrix also should contain one or more 
columns for interferents which are other analytes and stray 
features not of direct interest, and a column representing 
bacJcground. Vertical position in the column represents pixel; 
there generally is a separate submatrix for each subarray to form 
the Titatrix model, although one large matri« for the subarrarys is 
an alternative. The Ivaldi patent also teaches the inclusion of 
one or more mathematical derivatives in the matrix; however, as 
such derivatives are taught therein for correcting for spectral 
shift, they are not necessary in the present case which achieves 
the correction differently. Otherwise, however, the use of the 
matrix herein is essentially the same (but with different 
symbols) : 
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[C) - Cii](Uj 
where: 

M is a model matrix of known analyte data a^, background 
data bk and interferant data (where k is l to 
16) ; 

in R the superscript T indicates transverse operation, and 
the superscript -l indicates inverse matrix operation 
within the parentheses; 

U is a vector of spectral data u^ for a test sample; and 

G is a computed vector. 

This escample is for 16 pixels. The vector C includes a parameter 
cj representing a quantitative measure of a corresponding element 
in the test sample. The model is applied to the spectral data of 
the test sample using multiple linear regression by least 
squares . 

This method with a matrix model (without the derivative of the 
aforementioned Ivaldi patent, and without modification according 
to the present invention) is basically similar to Kalman 
filtering* for example as disclosed in U.S. patent Ho, 5,21B,5S3 
(de Loos-Vollebregt et al) , and in an article "Kalman Filtering 
for Data Reduction in Inductively Coupled Plasma Atomic Emission 
Spectrometry" by E.H. van Veen and M.T.C. de Loos-Vollebregt, 
Anal. Chem 63 1441 (1991) and previous articles referenced 
therein. Thus the present invention may be adapted to the Kalman 
filter format. 

According to a preferable aspect of the invention (FIG, 5) , 
before application of the matrix model to test spectral data, the 
matrix model 84 is shifted to the time (time 3) 80 of the 
acquisition of sample data 76. Thus the interpolated spectral 
shift 7 2 is applied 8 6 to the matrix model to effect a shifted 
matrix model 88. This shifted model is then applied so to the 
sample data 76 (taken at time 3) . The result is a display 92 of 
the composition 94 of the test sample. For another test sample, 
the intearpolation and shifting of the model is effected for the 
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new time of the new data acquisition. Spectral data for the 
drift standard is updated as necessary for the new interpolation. 

In developing the matrix model, inasmuch as the offset could 
drift during the (potentially lengthy) archive collection phase, 
there should be periodic updating of the current system " 
conditions as determined by repeated use of the drift standard in 
the manner set forth above. This provides a pixel offset number 
for the spectrum at each subarray, to which all further 
measurements on that subarray will be referenced to determine 
drift. Thus, each spectrum (spectral data) gathered at a 
subarray will have an offset associated with it. If spectra of 
several analytes are collected between drift standard 
measurements, the amount of offset is apportioned by 
interpolation (e»g. linearly) to each of them as a fraction of 
the time elapsed. The base spectral data for the model may b« 
formed by shifting each spectrum bade to a common initial time 
according the offsets. Alternatively and preferably, the offset 
is saved with the archive model data and applied later. 

The obtaining of a matrix model (PIG* 7, with some numerals the 
same as in FIG. 5 for similar steps or operations) is similar to 
the initial procedures of PIC, 5. This starts with acquiring 42 
initial spectral data 9 6 for a drift standard 4 4 at selected 
array positions 38 (FIG. 4) at an initial time (time i) 43. The 
Initial spectral data are compared 4 6 to a zero postion 48 
(explained above) to obtain initial offset data 98, and the 
initial offset data are used to obtain (compute) 54 an initial 
offset function 100 defining an offset for any array position. 
Preliminary spectral data 103 then is acquired for a )cnown 
concentration of the selected analyte 104 at an interim time 
(time 3) 80. Subsequent spectral data 10 6 for selected array 
positions is acquired 42 at a subsequent time (time 2) 60 for the 
drift standard 44, the subsequent spectral data are compared 4 6 
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to the zero postion 4S to obtain subsequent offset spectral data 
108, and the subsequent offset data are used to obtain Cconipute) 
54 a subsequent, offset function 110 defining an offset for any 
array position. The difference function €8 computed 66 between 
the initial offset function and the subsequent offset function is 
utilized by interpolation 70 to obtain an interim spectral shift 
112 for the interim tine 80. The interiio spectral shift is 
applied 74 to the preliminary spectral data 102 to effect the 
archive spectral data 114 which is stored for the base matrix 
model 120 associated with instrument conditions interpolated or 
extrapolated back to the initial time. 

As explained above, the interim spectral shift may be obtained 
alternatively via differences between the initial and subsequent 
offset spectral data- Also, it will be recognized that the 
foregoing steps of VXG. 7 for obtaining the model are achieved In 
the manner of FIG- S by preliminarily identifying the first time 
as the initial time, the second time as the subsequent tine, and 
the selected time as the interim time. An alternative in 
applying the shifts is to store all of the collected data for 
future computations of the shifts at the time of application of 
the matrix model to test sample data; the steps and means set 
forth herein and in the claims include this alternative. 

Upon query 115, if there are more analytes 116 (and there 
normally are many) , this procedure is repeated for each of these 
as well as interferents and the background until it is determined 
that the collection of archive data for the model matrix 12 0 is 
finished. At each drift measurement, if the drift is determined 
to be less than a preset standard, no correction is made. It is 
not necessary to use the drift standard after every analyte, and 
judgment should be exercised as to how often to use it, depending 
on early rate of drift or ether known factors on instrument 
stabilization, particularly with regard to temperature 
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stabilization. It is advantageous to collect and save each set 
of spectral data along with offset corrections interpolated baclc 
(in time) from -periodic data acquisitions with the drift 
standard . 

When the archive data •determinations are finished lie, a final 
acquisition nay be laade with the drift standard for the final 
interpolations/ and the computations are then made for the base 
matrix aodel 120. Alternatively, the conditions of the last 
spectral shifts may be used routinely without waiting for a 
processing from a later drift measurement, in such a situation, 
a later drift measurement may be made as an estimate of the 
degree of acceptability of the real-time results. In other 
words, a final result is computed and later validated upon 
measurement of drift to extablish whether such result has bean 
within tolerances. 

The spectral data are stored by pixel positions in the 
predetermined spectral Increments of pixel size. It may be 
advantageous to utilize pixel position (or pixel number) to 
replace units of wavelength (or wovenumber) , where the final 
output of the instrument is compositional information, with the 
spectral information being unnecessary. A small shift such as in 
drifting (e.g. from temperature change) will manifest as a small 
change in data levels in adjacent pixels. For the archive data 
collection of base spectral data for the matrix model, greater 
sensitivity is desired- For this, the instrument preferably is 
operated in a "high density" configuration, achievable by slit- 
scanning at a resolution of, for example, 4 intervals per pixel, 
whereby the entrance slit 17 (PIG. I) is positioned incrementally 
123 by a stepper motor 125. This enables the collection of 
reference spectra at a point spacing suitable for further 
mathematical manipulation. 
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With reference to PIG. when the model is applied 90 to test 
spectral data, the base spectral data 84 of the model is shifted 
88 as set forth, above. However, it is only necessary to use the 
model data in the pixel increments, e.g. every fourth data point 
in the higher resolution model. Thus, in suioraary, although the 
test spectral data are acquired at predetermined increments 
(pixels) . in associated subarrays, the preliminary spectral data 
for the model is obtained (and shift corrected for drifting 
during model archive acquisition) for a multiplicity of 
increments (e.g. 4) within each predetermined increment, and the 
shifted matrix modal is applied with model spectral data chosen 
at the predetermined increments. 

Initially, a large data base, or archive, of spectra encompassing 
all possible analytes, inter fcrents and background anticipated 
for the unknown test samples is collected. These spectra 
generally are obtained one at a time for pure element standards 
at concentrations where the maximum signal to noise ratio is 
obtained for the spectra, as well as maximum revelation of 
spectral features, but without overlaps from different analytes. 
Typically 60 to 70 analytes of individual atomic elements (or 
selected molecules for molecular spectroscopy) may be measured to 
assemble the archive for the matrix model. The individual 
analyte samples used for the archive are not necessarily the same 
as the drift sample which generally is a collection of several 
analytes as described above. 

All of the archival spectra could be adjusted to a common 
instrument condition. For example, the initial condition - which 
by definition has zero offset - might be chosen. Then, all 
spectra with finite offsets are shifted to the pixel they would 
have been centered on at the starting condition. Alternatively 
and preferably, the archival spectra and associated offset can be 
stored as-collected, and any shifting required can be deferred to 
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the offset (drift) conditions defined at run time. Critical to 
its success is the ability to accurately represent the spectra of 
all of the coioppnents in the mixtures to be analyzed. Since 
these must be in correct wavelength registry (or pixel location) , 
an optiiiiua conf igvuration occurs when all of these spectra are 
shifted to the conditions at which the run-time samples are 
collected, as determined by periodic use of the drift standard. 
Note that the run-time data does not need to (but may) be slit- 
scanned . 

The offset corrections with the interpolations are applied using 
a conventional method such as the Savitzky-Golay technique taught 
in an article "Smoothing and Differentiation of Data by 
Simplified Least squares procedures" by A. Savitzky and M.J.B. 
Golay, Analytical Chemistry 2&, 1627 (July 1964); five points are 
normally sufficient in the present example with the four interval 
slit scanning. After interpolation of the data the matrix model 
is assembled with every fourth point to correlate with pixel 
positions; i.e. with the first, fifth, ninth, etc, points. Other 
interpolation options include the conventional Lagrange and 
Fourier techniques. 

Alternative Drift Standard 

For the drift standard, an alternative to a standard sample of 
analyses is an optical interference element for effecting a 
series of regular spectral peaks 122 (FIG. 8), such as a fringe 
pattern with a multiplicity of interference fringes or secondary 
peaks spaced regularly across the spectral band including 
portions in each subarray 26. Use of a fringe standard in the 
form of a "low finesse etalon" is disclosed in the aforementioned 
U.S. patent No. 5,303,165 ("Ganz patent"), the portions relevent 
to such etalon being incorporated herein by reference. Such a 
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fringe standard is used in an Qmbodinent of the present 
invention, with modified procedures as explained below. 

Each secondary p^ak has an integer order number identifiable to a 
peak wavelength. A suitable fringe standard is a conventional 
Fabry-Perot low finesse etalon, advantageously an uncoated fused 
silica plate positioned ahead of the entrance slit to the 
dispersion system (e.g. gratings). ("Finesse" is defined 
conventionally as the ratio of spectral line spacing to line 
width in a series of regularly spaced interferoneter fringe 
lines.) In the optical train (FIG. 1) a source of broad band 
("white") light, e.g. from an incandescent filament or a sample- 
free plasma 18, is focused to the slit 17 by the input lens IS. 
For proper collimation, the etalon 124 (Pia. 9) in the optical 
train may be sandwiched between two proximate lenses 12«, 128, 
and this group is positioned between the input lens and the slit. 
Rays from the first, concave lens 126 are rendered parallel, 
passed through the etalon, and reconverged with a second, convex 
lens 126 to be focussed on the entrance slit. 

The low finesse etalon transmits the radiation in an interference 
fringe pattern forming a multiplicity of fringe peaks. The 
element has an index of refraction nCa,T> related to wavenumber a 
(reciprocal wavelength) and temperature T in a well-defined and 
wall known relationship which generally is ascertained from 
handbook tables or from standard equations such as those of 
Sellmeier discussed in "Temperature-Dependent Sellmeier 
Coefficients and Chromatic Dispersions for Some Optical Fiber 
Glasses", by G. Ghosh, M- Endo and T. Iwasaki, J. of Lightwave 
Technology, 12, 1338, August 1994. 

An integer order number m is identified to each peak wavenumber 
in accordance with the etalon thickness t by the standard fringe 
equation: 
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a=3ii/C2a(a,T) t] Eq. la. 

where the parentheses indicate dependence of the refractive index 
n on wavenumber a and temperature T. off-peak the order number 
is a non-integer. As an effective thickness t can be determined 
very accurately by procedures disclosed herein, the actual 
thickness need not be known accurately. Also, the procedures 
will average out nonuniform thickness and temperature-induced 
changes and account for slight misorientation of the eta ion in 
the instrument, which may vary from time to time. More broadly, 
following nonenclature of the Ganz patent. Eg. 1 is a correlation 
function and 2t is a correlation constant. 

The instrument is operated (PIG. lOA) with the low finesse ("LF") 
eta Ion t2A to acquire 131 secondary spectral data 129 
representing the fringe ("secondary") peaks which have peak 
locations 122 (FIG. B) ascertained with respect to pixels. Eg. 
la is utilized in conjunction with data also acquired 42 for an 
accurately known wavelength of a known peak 130 (or its centroid 
equivalent) such as a well known sharp absorption line from a 
known standard sample 133, a special lamp in the instrument, or 
an internal wavelength standard (e.g. Nd:VAC crystal taught in 
the Ganz patent) substituted for the low finesse etalon 124 in 
the optical train of FIG. 9 to produce primary spectral data 132 
representing the well defined primary spectral peak 130 having an 
Identified (absolute) wavenunber- The primary peak must be in 
one of the subarray positions which becomes a calibration 
subarray and conveniently is one of the selected subarray 
positions. The peak source operates as a calibration source to 
accurately locate the fringe peaks initially with respect to the 
pixel locations. 

Up to this point the procedure is sinilar to that of the 
aforementioned Ganc patent which is utilized with a continuous 
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array detector and discloses a ralatod procedure for determining 
etalon thickness, with fringe order number deterroined by counting 
fringes. Although described herein with respect to peaks, fringe 
valle/e may be used instead as disclosed In the Ganz patent and 
may be deemed peaks for the present purpose. If the present 
invention is applied to a continuous array detector, the further 
procedures of the Cans patent for determining order number, 
etalon thickness, and consequent peak wavelength may be followed. 
However, a fringe count cannot be achieved for a detactor with 
segmented subarrays for which the present invention is useful, so 
an alternative procedure is needed. 

In such alternative, present procedure, the nearest fringe peak 
134 that overlaps the pixel in a subarray containing the standard 
peak having wavenumber is identified 136 and fitted to 
determine the precise phase of the fringe peak relative to the 
standard peak. The location of this part of the fringe at the 
line peak is at an order number m|+«n where m, is an integer (for 
the near fringe peak) and *m is a fractional part of the order 
number. Eg. la becomes: 

B,+5m=:2n(aj,T)tao Eg. lb 

An initial estimate 138 of etalon thickness t, is made with a 
specified accuracy range *t for the possible thickness range, 
providing an initial minimum thickness t=t,-«t. An initial Sm 
also is estimated 140 (from the phase), allowing calculation 142 
from Eq. lb and the known a^, of a tentative (generally not an 
integer) for the minimum thickness. A nearest integral value is 
determined 14 4 for m|, and a corresponding thickness t then is 
calculated 146 from Eq. la, providing a first t,n| pair 147. 
Subsequent integers m, are listed 144 with corresponding 
thicknesses t calculated up to the maximum thickness t=t,4jt. 
For example, a^=3. 7707391x10* cm"*, n=l. 50023 for fused silica at 
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this wavenumber and tentperature T==26''C, thickness is estimated 
t.«500 iim with 5t-10 >4m, and initial 5m=0.25 fiti. Than ^or t-490, 
aii'«5543.586 with a nearest integer ino«5544 and a corresponding 
t-49 0.0366. The next integer is »,=5545 with a calculated 
t-490.125. This is continued vintil t=510 (approximately), 
resulting in a set of an initially large number of wavenunber- 
thiOcness (in„t) pairs 147. It will be recognized that the first 
calculation may be made for any thickness in the estimated range 
+/-5t^ and subsequent integers listed for thickness above and 
below the first thickness; with respect to the clains, such 
variation in procedure is to be deexied equivalent to that 
described above. Since the index of refraction changes faster in 
the deep UV spectrum, it is better to perforn the thickness 
noaeuremonts at higher wavelengths where temperature Induced 
variations from normal table values will be less significant. 

To narrow the choice (PIG. lOB) , another subarray is selected. 
For this, another line at another wavelength may be used, but 
only semi-accuracy is necessary, being obtainable from nominal 
instrument calibration or design, or the equivalent such as a 
standard sample of known analytes. Accuracy to one pixel is 
sufficient. In the new subarray, secondary (fringe) data ISO is 
acquired 148, a sinusoid is fitted to the this data, and a fringe 
peak 152 is selected 154 in a pixel having a wavenumber c known 
from the nominal calibration or a known peak in the subarray. 
Thickness values t from the first list of pairs are used to 
calculate 156 corresponding peak order numbers m within a 
fraction of a fringe uncertainty. An uncertainty of one pixel 
corresponds to a fraction of a fringe that can be calculated 
after the fit. The m values estimated from the thicknesses tried 
(that had passed from the previous stage) all must be an integer 
to within this fraction of a fringe to pass the current test. 
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since only integral order numbers ttj are correct, non- integers 
itfO are determined 158 to eliminate many pairs, and the remaining 
t^ra, pairs 164 are substantially reduced. (»<Non- integer" is to a 
prcdetemined level of precision.) This process is repeated 166 
(if necessary) for one or more other subarrays until only one 
pair is found 168, thereby establishing the effective etalon 
thickness t (or correlation constant 2t) • Corresponding order 
numbers ii| can then be selected and precise corresponding 
wavenunber^ calculated from Eq. la. It has been found that only 
three subarrays may be needed to converge to a unique answer. 
Once thickness is known, several known lines are measured to 
account for variations in the index and insufficient colligation. 

The true thickness of the etalon is determinable only to the 
e)Ctent that both the index of refraction and the calibration line 
are precisely Icnown. Even when both parameters are assumed to be 
accurate, uncertainties in the optical setup such as insufficient 
collination at the etalon may slightly alter results. A manner 
of correcting for such uncertainties, including the slight error 
possible in the thickness determination, is to use the thickness 
t determined as set rorth above, and an index of refraction that 
has been adjusted so that n,=i>.1-5n where n. is an adjusted index 
of refraction derived from the initially estimated index n, and a 
correction Sn to the index, and is dependent on wavenumber. A 
way to implement such a correction is to use a small number of 
additional standards with known emission lines and compute in 
values for each of these lines so that the order number will be 
what is actually observed for an estimated thickness. This 
correction may be plotted for a reasonable spread of lines, and a 
continuous function may be fitted to the data so that indices may 
be deducted for regions not directly measured. This correction 
is likely to be most important in the deep ultraviolet [UV) 
region where the index of refraction for silica increases 



(82) 



10-142054 



significantly as wavelength decreases, A mercury lamp is 
suitable for producing lines in the range of interest. 

With an etalon or the like substituted for a standard sample of 
selected analytes, the procedures for drift corrections are the 
sane as for the standard samples, both in assembling archive data 
for the basic matrix model and in shifting the model to current 
conditions. Selected fringe peaX data are identified as certain 
spectral data. Thus, in effecting the base spectral data for the 
model, initial spectral data are acquired at an initial time for 
the selected sinusoidal peaks with the etalon for the selected 
subarraye, the initial data are compared to the zero position, 
preliminary spectral data for tbe model is obtained for known 
concentrations of selected analytes, tgubsequent spectral data are 
acquired with the etalon at a subsequent time and compared to the 
zero position, and the preliminary data are shifted by 
interpolation with offset functions back to the initial time. 
Then, before and after data acquisition for test samples, first 
and second spectral data are taken with the etalon drift 
standard, spectral shift is determined, the matrix model data are 
shifted accordingly to the time of the test data, and the model 
is then applied to the test data- In recognition of the 
"absolute" nature of the etalon data, for this embodiment of the 
invention it may be desirable to reference the spectral data to 
wavenumber (or wavelength) rather than pixel position as for 
analyto drift standards. In this case, comparisons with zero 
position would be omitted. 

Temperature correction 

In addition to changes in the spectrometer itself arising from 
thermal expansion or contraction of the materials forming the 
instrument, temperature also affects the etalon by changing the 
index of refraction and thereby the fringe pattern- Two sets of 
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irfonaation are required, namaly am estinato of temperature at 
the etalon, and a degree of correction in the refraction index as 
a function of temperature. These may be determined in a series 
of steps. Data are taken at two levels of operational 
temperature (at least at the eta Ion ) , e.g. 2*C apart, held 
constant at each level. As the eta Ion will have changed 
thickness, its effective thickness should be determined at each 
teDperature by the foregoing procedures (taking into account that 
the previous value probably is close, o.g, within 1 ftm compared 
to 10 pm accuracy originally) . Despite the change in index of 
refraction, an estimated index n, may be used for both thickness 
determinations . 

At both temperatures the procedure is effected for n, by taking 
data for multiple lines, and two curves (or corresponding data in 
computer format) are derived to describe the effective index of 
refraction as a function of wavelength for both the initial and 
subsequent locations of the spectral peaks at the two 
temperatures. The peaks should include at least one line in the 
deep W where the index is more sensitive to temperature. This 
UV line serves as a "thermometer" since it shifts relative to the 
etaXon peaks with temperature, moving to an apparent wavelength 
away from its expected value. The initial and subsequent index 
curves are subtracted to provide the change in index as a 
function of wavelength for the temperature change. The changes 
(over wavelength) are divided by the degree of spectral shift for 
the peaks in the deep UV, to effect a set (or function) of slopes 
representing changes in index over wavelength relative to the 
shift in the V temperature" line. This information is stored for 
future use* Then, in ordinary operations, the "thermometer" line 
is observed periodically/ its shift is determined, the change in 
refractive index Sn is calculated from the slopes, and the 
changes are added to the nominal index n. to establish the 
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adjusted index of refraction n for the atalon as a function of 
wave length. 

A procodura <riG. 11) raay be utilized for temperature correction 
where the correlation function is dependent on index of 
refraction of the interference element, the index is dependent on 
temperature and subarray position, and the primary peak has a 
spectral position representative of temperature. While 
maintaining the interference element at a first temperature 170, 
the procedures described with respect to PIGS. XOA S lOB and Eq. 
la are effected 172 with a predetermined nominal index of 
refraction to ascertain first primary spectral data defining 
an associated first primary peak position (e.g. having wavenumber 

, a first effective correlation constant Ztj and associated 
order numbers »[- First values of the index of refraction n, are 
computed 174 with the correlation function <Eq. la) using the 
first effective correlation constant and associated order 
numbers, the first values being dependent on the nominal spectral 
positions for the selected subarrays- While maintaining the 
interference element at a second temperature, the procedures 
described with respect to riGS- lOA « lOB are repeated 172 with 
the nominal index to ascertain second primary spectral data 
defining an associated second primary peak position a,, a second 
effective correlation constant 2tj and associated order numbers 
n,- Second values of the index of refraction are computed 174 
with the correlation function using the second effective 
correlation constant and associated order numbers, the second 
values being dependent on the subarray positions. A value 
difference Dj-n, between each first value of the index and its 
corresponding second value is computed 176 for each corresponding 
subarray position, a position difference CTj-Ci between the first 
prinary peak spectral position and the second priwary peak 
spectral position is computed 178, and a difference ratio H=(n,- 
ai)/(aj-a,) of each value difference to the position difference is 
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coinputed 180. The resulting oet of diffaranco ratios as a 
functioji of subarray position is stored 182. 

At any selected time in which data are ta)cen with the eta Ion at a 
generally unknown temperature, subsequent primary spectral data 
are acquired 184 to define an associated subsequent prinary peak 
position a,, A subsequent difference aj-^j between the subsequent 
primary peak position and the first primary peak position is 
computed 186 as a measure of temperture change. A multiplication 
product 5n=(c3-<y,)*R of the subsequent difference and each 
corresponding difference ratio 182 is computed IBS to effect 
corrections in the index of refraction, and totals n^+dn of the 
changes and the nominal index of refraction are computed 190 and 
stored 191 to effect a temperature-corrected index of refraction 
n, for each of the selected subarray positions. The corrected 
index is then utilized in the correlation function for computing 
the spectral position (or wavenumber or wavelength) for the 
selected secondary peak in each selected subarray, in accordance 
with FIGS. lOA & lOB. 

Magnification 

As set forth above, slit scanning or the like is desirable to 
provide more points than allowed by the ordinary pixel size. 
This is used particularly for developing archive data for the 
matrix model, while the test data are acquired ordinarily in full 
pixel increments. More generally, the test spectral data, as 
well as the first and second spectral data, are acquired at 
predetermined intervals (e.g. pixels) in associated subarrays, 
the preliminary spectral data are obtained for a multiplicity of 
sub-increments smaller than the predetermined increments, and the 
shifted matrix model is applied with model spectral data culled 
for the predetermined increments. 
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In slit scanning, the entrance slit 123 (FIG. 1) is shifted 
laterally by a stepper motor 125 to move a spectral feature by 
selected sub- increments (fractions) of a pixel on a subarray, for 
example, 4 sub- increments per pixel- Essentially equivalent 
scanning may be achieved by other neans such as moving a lens or 
reflector in an optical train using a fixed slit. A full 
spectrum for all pixels is usually taken at each sub- increment, 
and the spectra are combined later to effect the actual sub- 
increment spectral data. The total number of sub-increments must 
be either spaced integrally or spaced with known intervals, 
requiring stepper motor action to be selected very accurately, 
i,e. precise knowledge of the amount of slit sub- increment is 
important. Any error results in an effective axis shift for the 
data. A corresponding correction factor for spectral position, 
called "magnification", should be made as close to one as 
practical, or eirrors develop in measurements of wavelength 
shifts. If the axis shift or magnification can be determined, it 
can be used with the scanning means (for example, stepper motor 
intervals) to correct the slit-scanning interval. The 
magnification (correction factor) is applied to the base spectral 
data to effect axis-corrected base spectral data in the base 
matrix model. 

If a system has, for example, precisely AO steps of the slit to 
move one pixel, then for four sub- increments per pixel, data 
collected at 10 step intervals would be precisely 0.2 5 pixel 
apart. However, systems generally may have a slightly different, 
predetermined, non-integral number of slit scans per pixel, e.g. 
40 + < where « is a fraction. The four scans separated by 10 
steps then would not cover one pixel, so that the net effect on 
the overall spectrum of combination of the 4 slit*-scanned spectra 
would be a nonuniformity of the spacing, since the spectra would 
be located at displacements of 0, 10, 20 and 3 0 steps in each 
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pixel* To coincide axactly with the next pixel would take an 
extra e fractional step. By determination of the exact number of 
steps to cover one piicel, a more accurate representation of the 
slit-scanned spectrum may be obtained. 

There are several vays to determine the magnification which is 
used to correct the scanning sub-increments to effect an integral 
total. In one embodiment (fig. 12), the slit is scanned 
preferably in smaller sub- increments or intervals than for normal 
slit scanning, for example 2 or 5 or 10 times smaller than 
normal. The number of slit steps should traverse at least one 
pixel from a nominal start position. For example, for a normal 
slit scanning in four sub-increments, an integral, nominal total 
number of 40 intervals or steps per pixel may be selected 192 as 
an approximation of the predetermined number. Based on this, a 
set of auxiliary total numbers of steps above and below the 
nominal total is selected, the set defining a series of 
multiplicities of sub-increments. The multiplicities (also 
denoted "index numbers") range smaller and larger in number than 
the nominal total, e.g. from 30 to 50. The scanning of each 
series in this example may be in increment steps of two pixels, 
i.e. 30, 32, 34, etc. Up to 50 sub-intervals per pixel. 

Slit scans 193 of the subarray are made for each series of index 
numbers on a selected pixel that should have a prominent feature 
(preferably a peak) from a sample, to acquire corresponding 
spectral data. Each scan (e.g. for 30, 32 etc.) produces a 
separate data series (first data 194a, second data 194b, third 
data 194o, etc.) across the pixel. Each original data series is 
pixel-axis shifted mathematically 19 by one pixel to effect 
first, second, third, etc. shifted data series 197. This data 
then is subtracted 198 from its corresponding shifted data (or 
equivalent vice versa) . Differences 199. may be presented 
actually as-is or, advantageously, are determined in a form of 
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standard root sum of squares (rss) of tha actual differences. As 
the index number approachee the actual number (e.g. approximately 
40) of scan steps to cover one pixel, the rss decreases. For 
example if the actual number is exactly 40, the rss would be zero 
for that index number and positive for 39 and 41* However, in 
general, there may not be an integral number of scan steps per 
pixel. 

To ascertain the actual (generally non-integer) number of scan 
steps, the rss data vs, index numbers are fitted 200 to two 
straight lines. One line is determined with index numbers below 
the nominal index (e.g. 40), and the other line is determined 
with index numbers above the nominal index. These two lines will 
have opposite slopes and will intersect 3 02 each other at the 
point of zero rss. The lines also should cross the axis at the 
same point but, if there is a small error, the nean location may 
be selected. The point defines the actual number 203 of slit 
scans per pixel. The ratio of the actual number to the nominal 
number defines the magnification. 

A check may be made with another pixel, advantageously the 
adjacent pixel at the nominal start position. The scans in this 
case thus may be made in the opposite direction beginning at the 
nominal start position. 

The foregoing, although effective, is slow due to the large 
number of scans. An alternative approach utilizes ordinary 
spectral data collection, with a derivative correction to 
wavelength shift. In a simple Taylor series approximation: 



y (x+5x) =y (X) +5x (dy/dx) 



£q. 2a 
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where y is spectral data and x is wavelangth (or equivalent pixel 
or sub- increment position) . In a discrete numerical embodiment 
for a spectrum. there is a function: 

Vi = aY, + l>(dVo/dXfl) + e Eq. 2b 

Where the Y's are spectral data vectors, is newly measured 
spectral data (possibly shifted by instrument drift) for a 
spectral feature, Yq is conveniently selected as the data for the 
original xinshifted archival data for that feature such as in the 
natrix model, and Xo is spectral position in a subarray. The 
tteasured data (subscript "1") conveniently (but not necessarily) 
are unscanned pixel data, whereas the archive data (subscript 
"0") may be in a higher (slit-scanned) density in sub- increments 
(e.g. 4 per pixel). (Data for could be data at pixel 
intervals of the archival data,) To apply the equation the 
archive data should be "thinned" to the same density as the 
measured data (one point per pixel) for data fitting to the 
equation. The derivative dY^/dXo »ay be taken at the high 
density and then thinned, or may be taken from the thinned data. 
The parameters a and b together represent the shift in spectral 
position (wavelength) as a parameter ratio b/a. The parameter e 
is the fitting residual. 

The axis shift associated with an axis magnification deviating 
from one may be determined by trial and error. A tentative 
magnification is estimated, advantageously assumed to be 1.00. 
The apparent shift b/a is computed with Eq. 2b for a plurality of 
subarrays, for example the same selected subarrays used for 
determining the offset data so that no new data are required. 
When the correct magnification is used, these shifts are 
proportional to the size of the slit^scan steps. With a straight 
line fitted between b/a and the size, for all of the selected 
subarrays, the root-raean-square (rms) of the fit gives an 
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indication of the deviation from the correct magnification. This 
ras raay be large at the beginning before the correct 
magnification is approached- 

To correct this a corrected magnification is estimated, e.g. 
varied from a. 00. The size of spectral position Xo sub- 
increments are adjusted by the magnification to effect adjusted 
spectral data. Eq. 2b is applied again, and the parameter ratios 
D/a vs. step size are refitted to a straight line. This is 
repeated until the rms is minimized. Trial and error may be 
replaced by an algorithmic method such as the conventional 
biaectional search. Alternatively a number of calculations could 
be made with a range of magnifications and the one with the 
smallest error selected. The magnification generally will be a 
function of the n, y coordinates of the subarrays, as for the 
offsets. 

More specifically (FIG. 13), a tentative eize of the sub- 
increments (steps) 204, and a corresponding tentative totals 206 
of sub- increments in a pixel increment, are estimated 208 for 
each of the selected subarrays. These generally will depend on 
subarray on the x,y positions. Spectral data for the derivative 
is acquired 212 either at this time or, preferably, as previous 
archive data 114, and the derivative 210 is calculated 212. 
Current spectral data Y, is acquired 214 for each of the selected 
subarrays and tentative steps. The current data and derivative 
are fitted 216 to Eq. 2b to compute the apparent shift b/a. The 
latter vs. the tentative sizes 204 of sub-increments for the 
selected subarrays are fitted 218 to a curve to establish a 
deviation from a straight line, the "curve" advantageously being 
a straight line with computation of the root mean square (rms) 
220 of the deviations. Based on this rms, a corrected 
magnification 222 (variation from 1.00) is estimated 224. Ttiis 
is multiplied 225 by the previously estimated sizes 204 of sub- 



increments to yield new tentative sizes and totals 206. The 
sequence is repeated 226 until, on query 228, the rms is 
minimized to within a preselected level, at which time the last 
tentative total 208 is detemined 230 for the final total 232. 

while the invention has been described above in detail with 
reference to specific embodiments, various changes and 
modifications which fall within the spirit of the invention and 
scope of the appended claims will become apparent to those 
skilled in this art. Therefore, the invention is intended only 
to be limited by the appended claims or their equivalents. 
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4. Brief Explanation or tne ijravvlna^ 

FIG. 1 is a schematic drawing of a spectroraetric instrunent used 
for the invention. 

riG. 2 is a schematic view of the front surface of a detector in 
the instrument of FIG. 1, showing subarrays of photosensitive 
pixels. 

FIG. 3 shows the shape of a spectral feature and associated 
levels of photo energy received by pixels of the detector of FIG. 
2. 

FIG, 4 shows a subset of selected subarrays having position 
coordinates on the detector of FIG. 2. 

FIG. 5 is a flow diagram of steps and means for an enbodinent of 
computations associated with the instrument of FIG. 1 according 
to the invention. 

PIG, 6 shows a spectral feature associated with the flow diagram 
of FIG. S. 

FIG. 7 is a flow diagram of steps and means for another 
embodiment of computations associated with the instrument of FIG. 
1 according to the invention. 

FIG. 8 shows a series of regular spectral peaks from an drift 
standard utilized for an aspect of the invention, and a 
superimposed spectral feature, associated with the flow charts of 
FIGS. 5 and 7. 

FIG. 9 is a schematic diagram of an optical train incorporating 
the drift standard of FIG. 8 into the instrument of FIG. 1. 

FIGS. aOA and lOB are flow diagrams of steps and naanc for 
utilizing the drift standard of FIG, 8 according to the 
invention. 

FIG. 11 is a flow diagram for temperature correction utilized in 
conjunction with FIGS lOA and lOB, 

FIG, 12 is a flow diagram for adjusting magnification in a slit 
scanning embodiment incorporated into FIG. 1. 

FIG. 13 is a flow diagram for adjusting magnification in an 
alternative slit scanning embodiment incorporated into FIG, 1. 
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ABST3UVCT OF THE DISCLOSURE 



A spectr ©metric instrument includes a detector with detecting 
subarrays on small portions of the surface. Spectral data ore 
acqtuired for selected subarrays at a first time for a drift 
standard, and compared to a zero position to obtain first offset 
data. Data are acquired similarly at a second time to obtain 
second offset data. The offset data are utilized to obtain a 
spectral shift for any subarray position at any selected tirte. 
The shift is applied to a matrix model used for converting test 
data to compositional information. 



